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Preface 


The  purpose  of  this  study  was  to  determine  if  several 
proposed  pilot  induced  oscillation  (PIO)  prediction 
techniques  could  predict  longitudinal  PIO  tendency  in  a 
variety  of  aircraft  configurations  prior  to  flight  test.  The 
techniques  were  first  applied  to  an  existing  data  base  to 
provide  insights  into  their  application  and  limitations. 

Based  on  the  analytical  results,  18  aircraft/flight  control 
systems  were  selected  for  flight  testing.  Each  technique  was 
used  to  predict  the  PIO  tendencies  of  the  configurations 
prior  to  flight. 

Although  the  flight  test  results  were  inconclusive  due 
to  the  limited  number  of  configurations,  the  results  provided 
some  intriguing  possibilities  for  future  PIO  research.  The 
data  base  established  here  is  excellent  for  use  in  future 
research;  the  experiment  was  tightly  controlled  and  provided 
extremely  consistent  results. 

The  joint  AFIT/USAFTPS  program  under  which  this  research 
was  conducted  provided  a  unique  opportunity  to  apply  academic 
research  to  an  operational  test.  In  performing  both  my 
research  and  the  flight  test  program  I  was  aided  by  many 
people.  I  must  especially  thank  my  thesis  advisor,  Dr.  R.  A. 
Calico,  for  providing  me  with  the  topic  in  the  first  place. 
Also,  much  thanks  must  be  given  to  Major  J.  T.  Silverthorn, 
my  advisor  at  TPS,  for  providing  invaluable  assistance  during 
the  flight  test  program,  particularly  the  data  reduction.  I 
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am  also  deeply  idebted  to  Ralph  Snith  (now  with  NASA)  for 
reviewing  much  of  my  original  research  and  providing  me  with 
seme  excellent  insights  to  my  approach.  Much  credit  must 
also  be  given  to  my  test  team  at  TPS.  The  members  included 
Capt  Dave  Eidsaune,  Capt  Rick  Bennett,  Capt  (Japan)  Seiichi 
Miyamoto,  and  lLt  Carter  Wilkinson.  By  pulling  together,  we 
somehow  managed  to  complete  a  six  week  test  program  in  a 
little  over  two  weeks.  Likewise,  I  must  thank  all  of  the 
Calspan  safety  pilots  involved  in  the  program  (Mike  Parrag, 
Bob  Harper,  and  Russ  Easter)  and  the  Calspan  maintenance 
crews,  who  worked  many  long  hours  and  weekends  to  allow  us  to 
complete  the  test  program  on  schedule.  I  owe  a  special 
thanks  to  Mike  Parrag  and  Lou  Knotts  (also  of  Calpsan)  for 
their  constructive  comments  and  assistance  during  the 
planning  phase  of  the  flight  test  program.  Their  assistance 
undoubtedly  helped  us  avoid  many  of  the  pitfalls  which  often 
plague  research  programs  of  this  nature.  Finally,  I  must 
thank  Lt  (USN)  Ron  "Weasel"  Weisbrook  for  taking  the  time 
from  his  busy  schedule  to  proofread  several  of  my  drafts. 


Eileen  A.  Bjorkman 
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""-The  purpose  of  this  study  was  to  determine  if  pilot 
induced  oscillations  (PIOs)  can  be  predicted  prior  to  flight 
using  existing  PIO  prediction  techniques.  Two  techniques  to 
predict  longitudinal  PIO  tendencies  (Ralph  Smith's  theory  and 
Roger  Hoh's  bandwidth  method)  were  studied  analytically  using 
an  existing  PIO  data  base.  Suggestions  were  made  for  both 
techniques  to  allow  prediction  of  PIO  rating.  The  two 
techniques  were  then  applied  to  18  aircraft/ flight  control 
system  landing  configurations.  The  18  configurations  were 
flight  tested  using  a  flared  landing  task  with  the 
USAF/Calspan  variable  stability  NT-33A.  Smith's  theory 
correctly  predicted  the  PIO  tendencies  and  frequencies 
provided  the  configuration  was  not  sensitive  to  the  pilot 
model  used.  A  suggested  modification  to  Snith's  theory 
correctly  predicted  PIO  ratings  within  an  average  of  0.6 
rating.  A  suggested  modification  to  Hoh's  bandwidth  method 
predicted  PIO  ratings  within  an  average  of  0.5  rating. 

The  limited  data  base  was  too  small  to  draw  any  definite 
conclusions.  Recommendations  for  further  study  included 
collecting  more  PIO  data  and  using  existing  data  bases  and 
simulator  studies  to  better  define  the  two  techniques  and  to 
gain  physical  insights  into  PIO  mechanization,  i 
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FLIGHT  TEST  EVALUATION  OF  TECHNIQUES  TO 
PREDICT  LONGITUDINAL  PILOT  INDUCED  OSCILLATIONS 


I.  Introduction 


Background 

A  pilot-induced  oscillation  (PIO)  is  an  inadvertent, 
sustained  oscillation  of  a  pilot/vehicle  system  (1:1). 
Characteristically,  a  PIO  is  difficult  or  impossible  for  the 
pilot  to  stop  unless  he  removes  himself  from  the  loop,  i.e., 
the  aircraft  is  stable  both  stick-fixed  and  stick  free  (2:2). 
PIO's  typically  occur  during  tasks  for  which  the  pilot  is 
particularly  concerned  about  tight  control  of  the  aircraft, 
such  as  during  landing,  takeoff,  refueling,  and  formation 
flying . 

PIO's  have  been  documented  since  the  beginning  of  manned 
flight.  Although  even  the  Wright  Flyer  had  a  mild 
longitudinal  oscillation,  PIO  did  not  become  a  serious 
problem  until  high  performance  jets  emerged  in  the  fifties. 
The  high  speeds  and  fully  powered  control  systems  of  modern 
aircraft  can  be  potentially  lethal?  relatively  small  pilot 
inputs  can  cause  a  rapid  buildup  to  catastrophic  loads  (1:1). 
A  well  documented  PIO  occurred  during  the  early  sixties  when 
an  early  version  of  the  T-38A  sustained  a  severe  PIO  (seven 
cycles  of  ±4  g's  building  up  to  ±8  g’s)  after  the  pilot  shut 
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down  a  malfunctioning  pitch  damper  with  the  elevator 
mistrimmed  (1:1).  This  PIO  is  clearly  demonstrated  by  the 
time  histories  in  Figure  1.  The  problem  was  traced  to  the 
bobweight  in  the  flight  control  system,  which  was 
subsequently  modified  to  prevent  a  reoccurrence. 

Typically,  PIO's  have  gone  undetected  until  the  final 
flight  testing  or  early  production  stages  of  an  aircraft  when 
aggressive  pilot  behavior  is  more  likely.  Since  fixes  to  an 
aircraft  at  this  late  stage  are  difficult  and  expensive, 
there  have  been  many  attempts  to  analytically  predict  PIO 
tendencies  in  an  aircraft/flight  control  system  before  it  has 
been  built.  Unfortunately,  the  current  military  flying 
qualities  specification,  MIL-F-8785C  (3),  provides  no 
guidance  for  precluding  PIO  by  design.  Although  two  possible 
PIO  prediction  techniques,  the  R.  smith  criterion  (2)  and 
Hoh's  bandwidth  method  (4),  were  proposed  for  inclusion  in 
MIL-F-8785C,  the  current  version  simply  states  that  an 
aircraft  will  not  PIO. 

Objectives 

The  purpose  of  this  study  is  to  evaluate  and  attempt  to 
refine  some  of  the  more  recent  longitudinal  PIO  prediction 
theories.  The  ultimate  goal  is  to  use  these  theories  to 
predict  PIO  tendency  for  a  variety  of  aircraft/flight  control 
system  configurations  in  the  approach/ landing  phase  of  flight 
and  to  verify  these  predictions  using  the  Calspan/USAFTPS 
variable  stability  NT-33.  The  theories  to  be  examined  are  R. 
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TIME  -  SECONDS 


anith's  longitudinal  PIO  criteria  (2,5)  and  R.  Hoh's 
bandwidth  method  (4).  Although  Hoh’s  method  has  been 
proposed  as  a  handling  qualities  predictor,  it  will  be 
employed  in  this  study  to  see  if  it  can  be  used  to  predict 
PIO  tendencies . 

Specific  objectives  include: 

1.  Determine  if  any  parameters  obtainable  from 
aircraft/flight  control  system  dynamics  can  be  used  to 
predict  PIO  ratings  and  frequencies. 

2.  Determine  if  existing  PIO  prediction  techniques  need 
refinement  and  make  suggestions  for  improvement. 

3.  Use  each  technique  to  predict  prior  to  flight  the 
PIO  tendencies  of  a  variety  of  aircraft/ flight  control  system 
configurations;  determine  the  percentage  of  correct  PIO 
predictions  made  by  each  technique,  including  PIO  rating  and 
frequency. 

Approach 

Specifically,  the  approach  taken  will  be: 

1.  Study  the  two  techniques  and  develop  or  obtain 
computer  programs  to  ease  their  implementation. 

2.  Use  the  1978  Calspan  Landing  Approach  Higher  Order 
System  (LAHOS)  data  (6)  to  see  how  well  the  theories  predict 
PIO's  noted  during  this  study.  Determine  if  improvements 
need  to  be  made  to  any  of  the  techniques. 
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3.  Determine  if  any  parameters  obtained  from 
aircraft/flight  control  system  dynamics  can  predict  actual 
PIO  ratings.  Also,  determine  if  any  of  the  techniques  can 
predict  actual  PIO  frequency. 

4.  Determine  a  set  of  operationally  realistic  aircraft/ 
flight  control  system  configurations  for  the  approach/ landing 
phase  which  can  be  simulated  on  the  NT-33.  Configurations 
will  be  selected  to  verify  and  refine  the  techniques  based  on 
insights  obtained  from  examining  the  LAHOS  data. 

5.  Use  the  theories  to  predict  which  configurations 
from  step  4  will  be  PIO  prone  and  which  will  not  be  PIO 
prone.  Predict  PIO  ratings  and  frequencies,  if  possible. 

6.  Flight  test  the  configurations  using  the  NT-33. 
Confirm  or  refute  the  predictions. 

Scope 

It  should  be  noted  that  this  study  is  not  an  attempt  to 
develop  a  new  PIO  theory.  Rather,  it  will  attempt  to  verify 
and  refine  existing  theories  and  to  determine  whether  any  of 
the  theories  could  be  useful  for  inclusion  in  a  future 
version  of  MIL-F-8785C.  The  study  will  examine  only 
longitudinal  PIO  and  will  not  attempt  to  study  nonlinear 
control/feel  system  dynamics.  Also,  it  is  important  to  point 
out  that  the  two  techniques  used  in  this  study  are  not  the 
only  methods  currently  available  which  may  be  useful  to 
predict  PIO  tendencies.  These  two  theories  were  selected 
because  they  have  both  been  studied  previously  for  inclusion 
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in  versions  of  MIL-F-8785C.  In  addition,  the 
approach/landing  phase  of  flight  was  selected  for  analysis 
because  it  is  a  repeatable,  high  gain  task  which  provides 
more  consistent  PIO  data  than  "up  and  away"  flight.  Both  the 
theories  examined  are  directly  applicable  to  "up  and  away" 
flight. 

More  importantly,  the  intent  of  this  study  is  not  to 
deal  with  PIO  tendency  from  only  an  academic  point  of  view, 
but  to  deal  with  it  from  an  operational  perspective.  A  mild 
PIO  tendency  that  can  be  accurately  predicted  may  be  of 
academic  interest  only.  There  is  little  value  in  predicting 
a  mild  PIO  if  the  pilot  does  not  find  it  objectionable  and  it 
does  not  impair  his  ability  to  perform  a  task.  It  is  more 
important  to  find  a  theory  which  not  only  accurately  predicts 
a  PIO  tendency,  but  can  also  tell  us  something  about  the 
impact  the  resulting  PIO  tendency  will  have  on  the  aircraft's 
mission. 
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II.  Theoretical  Development 


This  section  first  presents  the  aircraft  equations  of 
motion  needed  for  longitudinal  PIO  analysis,  and  then 
summarizes  the  R.  Smith  and  bandwidth  theories.  Finally,  an 
example  will  be  given  which  applies  each  theory  to  the  YF-17 
as  simulated  using  the  NT-33  during  the  LAHOS  study. 

Equations  of  Motion 

The  equations  of  motion  developed  in  this  section  are 
for  the  vehicle  body  axes  using  "lumped"  stability 
derivatives  as  described  in  the  parameter  identification 
technique  of  the  LAHOS  study  (6:211).  The  "lumped"  stability 
derivatives  include  the  and  terms  in  the  appropriate 
stability  derivatives  (see  Appendix  A  for  further 
clarification) .  The  general  longitudinal  equations  of  motion 
are  presented  first;  the  specific  open  loop  transfer 
functions  necessary  for  longitudinal  PIO  analysis  are  then 
developed . 

The  Laplace  transformed  longitudinal  perturbed  equations 
of  motion  for  the  airframe,  referenced  to  the  vehicle  body 
axis  system,  using  "lumped"  stability  derivatives,  and 
considering  only  vertical  gusts  and  elevator  deflection  as 
inputs  are  (7:256) 
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[1] 


su  +  Woq  +  gcoseo0  -  Xuu  +  Xw(w-wg)  +  Xqq  +  Xfi 

sw  -  u  q  +  gsine  o  =  Z  u  +  z  (w-w  ) 
c>  3  o  u  w  g 

+  Zq(q+SWg/Uo)  + 

sq  *  M  u  +  M  (w-w  )  +  M  (q+sw  /U  )  +  M  £ 

^  u  w'  g'  q  ^  g'  o  6e  e 

where 

u  =  perturbed  forward  speed 
w  =  perturbed  downward  speed 
e  =  perturbed  pitch  angle 
q  =  perturbed  pitch  rate 
$e  =  elevator  deflection  input 
wg  =  vertical  wind  gust  input 
UQ  =  equilibrium  forward  speed 
WD  =  equilibrium  downward  speed 
0Q  =  equilibrium  pitch  angle 

Using  the  two  relations  q=se  and  a»tan_1 (w/UQ) *w/UD 
(valid  for  small  angles)  equation  [1]  becomes 

su  +  W  se  +  gcose  0  =  X  u  +  X  (U  a-w  )  +  X  so  +  X,  6e 
o  o  u  wo  g  q  oe 

saUo  ‘  V9  +  9sinV  =  Zuu  +  Zw(Uoa"V 

+  Z  ( S0+SW  /U  )  +  Z  6 e  [2] 

q  g  u  6  e 

s20  =  M  u  +  M  (U  a-w  ) 
u  w  o  g 

+  M  (s0+sw  /U  )  +  M  „  6e 
q  g '  o  6e  e 

where  a  =  perturbed  angle  of  attack. 

Usually,  X  =Z  =0  is  a  good  approximation  (7:273-277). 
Rearranging  terms  and  putting  [2]  into  matrix  form  yields 
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In  addition,  the  following  transformed  kinematic 
relationships  define  the  normal  acceleration  at  the  pilot 
station: 

se  »  q 

a  =  w  -  U  s  0  [4] 

z  o 

a  =  a  -  1  s^e 

Zp  Z  X 

where 

a  =  downward  (normal)  acceleration  of  the  aircraft 
z 

center  of  gravity 

a  =  pilot-felt  normal  acceleration  (normal  acceleration 

Zp 

at  the  pilot  station) 

1  =  pilot's  location  forward  of  the  c.g. 


Selecting  one  input  (set  either  6e  or  Wg=0)  and  using 
Cramer's  Rule,  the  following  "generic"  transfer  function  can 
be  developed: 
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Output 

Input 


Numerator 


Numerator 


(s-X  ) 
u 

-X  u 
w  0 

-Z  /U 
u'  0 

(s-Z  ) 
w 

-M 

u 

-M  U 
w  0 

(W  s+qcos8  ) 
o  o 


A  =  (s2+2c  «  s+u  2) (s2+2t__u  s+u  2) 
p  p  p  sp  sp  sp 


(S2-Mq3) 

The  characteristic  equation  of  the  system  (A)  is  a 
fourth  order  polynomial  which  factors  into  two  second  order 
polynomials  and  can  be  written  as 

'i 

[5] 

These  two  modes  are  usually  oscillatory  and  represent  the 
free  longitudinal  motions  of  the  aircraft.  They  are  called 
the  "short  period"  and  the  "phugoid."  The  short  period  mode 
is  a  relatively  well  damped,  high  frequency  oscillation, 
while  the  phugoid  mode  is  a  lightly  damped,  relatively  low 
frequency  oscillation.  The  short  period  mode  is 
characterized  by  small  changes  in  forward  velocity,  u,  and 
large  amplitudes  of  0  and  a. 

A  two  degree  of  freedom  short  period  approximation  to 
the  longitudinal  equations  of  motion  can  be  obtained  by 
assuming  constant  airspeed  (setting  u  to  zero).  Both  the 
fourth  order  dynamics  and  the  short  period  approximation  were 
initially  used  in  this  study  and  compared;  however,  for  the 
analysis  of  landing  configurations,  the  short  period 
approximation  in  general  did  not  give  satisfactory  results. 
The  poor  results  achieved  with  the  short  period  approximation 
may  be  due  to  the  constant  airspeed  assumption.  DiDomenico 
(8:16)  retained  the  phugoid  mode  for  his  landing  flare 
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handling  qualities  study  since  the  phugoid  accounts  for  the 
airspeed  "bleedoff"  during  the  flare.  In  the  following 
development,  the  complete  dynamics  of  equation  [1]  are  used. 
The  fourth  order  longitudinal  characteristic  equation  is 

t  m  S4  +  As3  +  Bs2  +  Cs  +  D  [6] 

where 


A  =  -M  -Z  -Z 
q  w  u 

B  =  X  (M  +Z  )+M  Z  -U  M  -X  Z  +W  M 
u  q  w  qw  ow  wu  ou 

C  -  -X  (Z  M  -U  M  )+Z  (X  M  +W  M  )-M  (U  X  +W  Z  -gcose  ) 
u  w  q  o  w'  u'  w  q  o  w'  u'  o  w  o  w  *  o 

+gM  sine 
3  w  o 


D  "  gcose  (Z  M  -M  Z  )-gX  M  sine 
ouwuw  uw  o 


The  numerators  of  the  transfer  functions  important  to 
PIO  study  will  now  be  developed. 


Pitch  to  Elevator  Deflection  Numerator. 


where 


N*j  (s)  = 
6e 


(s-X  ) 

-X  U 

u 

W  0 

-Z  /u 
u'  o 

(s-zw) 

-M 

-M  U 

u 

W  0 

N !  (s) 
fie 


+  B  s  +  C 
6  6 


[7] 


A  =  M 
e 


6, 


B  =  X .  M  +ZX  M  -M  (X  +Z  ) 
e  6e  u  6e  w  6  e  u  w 

C„  =  X,  (Z  M  -M  Z  )+Z,  (X  M  -X  M  )+M,  ( Z  X  -X  Z  ) 
0  6e  uwuw  6ewuuw  6ewiiwu 
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Angle  of  Attack  to  Elevator  Deflection  Numerator.  This 


transfer  function  is  needed  to  find  the  normal  acceleration 
transfer  function. 


(S-Xu) 

Xfi 

6  e 

(W  s+gcose  ) 

0  o 

^e1’1  * 

-ZJVo 

Z.  /Uo 
o  e  o 

(-Uos+gsin0o)/Uo 

-M 

u 

(s2-Mqs) 

N°  (s)  =  A  s3  +  B  s2  +  C  s  +  D 


where 


A  =  Z.  /U 
a  6  e  o 

B  =  [X,  Z  -Z.  (M  +X  )+M.  U  ]/U 
a  Je  u  6 e  q  u  oe  O  O 

Ca  =  [X,  (-Z  M  +U  H  )+Z.  (X  M  +W  M  ) 

»e  u  q  o  u  Qg  u  q  o  u 

+V(-9sinVXuVWoV]/Uo 

Dq  -  [gcos8o(MuZ5^-ZuMfs)+gsineo(M(eXu-MuX4(i)]/U0 


Pilot-Felt  Normal  Acceleration  to  Elevator  Deflection 
Numerator.  Fran  equation  [4] 


where 


2 

a  =  sw  -  sU  0  -  1  s  e 
zp  ox 

=  sU  a  -  sU  6  -  s2l  0 

O  O  X 

N^Zp(s)  =  sU  N®  (s)  -  sU  N®  (s)  -  s2l  N ®  (s)  [9] 

o  6e  06e  X  6  e 

N^zp(s)  =  A  s^  +  B  s^+C  s2  +  D  s 

6a  a  v  a  a  7  a? 

“  ^  7  _  ^  n  ^  n 
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% 

o  a 

x  e 

B 

=  U  (B 

-  A)  - 

1  B 

3ZP 

O  a 

0 

X 

C 

-  U  (C 

-  B.)  - 

1  c 

azP 

D 

azP 

0  0 

=  U  (D 

0  0 

0 

-  C0> 

X 

Pitch  Loop  Equivalent  Command  Gust  Numerator. 


Nw  (s)  = 
wg 


(s-X  ) 
u 

-X  u 

-X 

W  0 

w 

-VUo 

<s-Zw> 

-z  /u 

w'  o 

-Mu 

-Vw 

-M+Ms/U 
w  q  o 

where 


«£(»>-  *e„.*3  +  a2  +  ce„  s 


Jwt 


JW, 


Wr 


6w 


g 


M  /U 

q  o 


o 

Cn  =  X  (M  +Z  M  /U  )-X  (M  +Z  M  /U  ) 
0W  u  w  w  q'  o  w  u  u  q'  o 

O 


[10] 


The  primary  transfer  functions  needed  for  evaluating  PIO 

tendencies  are  the  equivalent  gust  command  transfer  function, 

e/w  ,  the  pitch  to  stick  force  dynamics,  e/F  ,  pilot-felt 
9  8 

normal  acceleration  to  stick  force  dynamics,  a  /F  ,  and 

z  p  s 

pilot-felt  normal  acceleration  to  pitch  rate,  a  /e.  The 

Zp 

last  three  can  be  derived  by  knowing  the  e/ 6  and  a  /d 

e  z  p  e 

transfer  functions  and  the  control,  feel,  and  actuator 
systems  of  the  aircraft  as  shown  in  Figure  2. 
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Figure  2.  Aircraft/Control  System  Block  Diagram 


There  is,  obviously,  no  "generic"  way  of  deriving  the 
feel  system,  control  system,  and  actuator  dynamics  (fi_/F_). 
Once  this  transfer  function  is  known,  however,  the  three 
necessary  transfer  functions  can  be  derived: 

e/F  (s)  =  [0/ 6(s)][</F(s)] 
s  e  e  s 

aZp/Fs(s)  =  [azp/<e(s,][ie/Fs(s)]  [12] 

a  / e ( s )  =  [l/s][a„  /«  (s)J[«  /e(s)] 
z  p  z  p  e  e 

Having  developed  the  transfer  functions  needed  for 
longitudinal  PIO  analysis,  the  theoretical  development  will 
turn  to  a  discussion  of  the  two  PIO  theories  used  in  this 
study. 
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Stoith's  Theory 


Ralph  Smith  identifies  two  types  of  PIO  in  Reference  2, 
and  a  third  PIO  type  in  Reference  5.  Type  I  PIO  is  initiated 
by  pitch  attitude  control  and  is  postulated  to  occur  in 
situations  where  the  pilot  "switches"  control  from  tracking 
pitch  attitude  to  tracking  pilot-felt  normal  acceleration. 
Type  II  PIO  is  initiated  by  abrupt  turbulence  or  nontracking 
abrupt  maneuvering,  such  as  SAS/CAS  start-up  or  shut  down  or 
trim  malfunction.  A  third  type  of  PIO,  which  will  be  called 
Type  III  PIO  in  this  report,  was  discussed  by  Smith  in 
Reference  5.  This  type  of  PIO  is  initiated  by  pitch  attitude 
tracking  only;  i.e.,  the  pilot's  acceleration  channel 
dynamics  are  irrelevant.  To  identify  an  aircraft  as  being 
PIO  prone  or  free,  all  three  types  of  PIO  must  be  examined. 

The  discussion  of  Snith's  PIO  criteria  presented  below 
is  taken  from  References  2  and  5.  For  each  type  of  PIO,  the 
theory  is  first  presented,  followed  by  a  summary  of  the 
assessment  rules.  At  the  end  of  the  section,  the  results  of 
applying  Smith's  criteria  to  the  YF-17  as  simulated  during 
the  LAHOS  study  are  presented. 

Type  I  PIO.  Smith's  technique  may  be  summarized  as 
follows:  There  can  be  a  frequency  at  which  the  power  spectral 
density  of  the  pilot's  normal  acceleration  due  to  pitch 
attitude  tracking  is  "sufficiently"  narrowband.  If  such  a 
frequency  exists,  there  is  a  high  probability  during  a  high 
gain  tracking  task  the  pilot  will  switch  from  tracking  pitch 
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to  tracking  normal  acceleration.  In  this  case,  the  frequency 

is  said  to  be  "subjectively  predictable."  A  suggested 

threshold  for  the  magnitude  of  normal  acceleration  the  pilot 

must  sense  to  attempt  acceleration  tracking  is  laZp/MwR)|  > 

0.012  g/deg/sec  (2:36).  Then,  if  the  phase  margin  of  the 

pilot-felt  normal  acceleration  to  stick  force  dynamics  is 

less  than  zero  at  the  subjectively  predictable  frequency,  the 

aircraft  will  have  a  tendency  to  PIO  at  that  frequency. 

To  understand  how  Type  I  PIO  can  occur,  Smith's  pilot 

model  must  be  presented.  The  model  is  shown  in  Figure  3. 

The  pilot  compensation  dynamics,  Yp(jw),  are  in  the  form  of 

the  servo  model  presented  in  Reference  9.  The  switch  in  the 

diagram  is  either  on  or  off;  i.e.,  the  pilot  either  tracks 

pitch  angle  or  normal  acceleration.  Under  normal 

circumstances  the  pilot  will  track  pitch,  but  Smith  proposes 

that  if  at  seme  point  the  power  spectral  density  of  a  due 

ZP 

only  to  closed  loop  control  of  pitch  attitude  (denoted  by 
*  )  is  sufficiently  narrowband,  the  pilot  will  switch  to 

3  Za  Z 

tracking  normal  acceleration.  In  order  for  *  to  be 

a  zaz 

narrowband  the  closed  loop  pilot-vehicle  system  for  pitch 
tracking  must  be  resonant,  as  wi!  1  be  discussed  later.  One 
of  Smith's  necessary  conditions  for  PIO  is  that  the  pilot  at 
some  point  begins  to  track  normal  acceleration.  The  theory 
describing  this  "switch"  is  beyond  the  scope  of  this  study, 
but  a  short  discussion  of  the  computation  is  presented  below. 
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K  (Tt jw+1) 


*  V*  '  T 

ij.)  =-£-!; 


,-’e3* 


(Tjju+l) 


K  (a  )  »  0  when  the  PSD  of  a  (t)  is  broadband 
a  z  p  Zp 

K  (a  )  =  K  when  the  PSD  of  a„  (t)  is  narrowband 
a  z  p  ao  z  p 

K  (a  )e-Ta3“  =  pilot's  acceleration  channel  dynamics 
a  z  n 


Figure  3.  Smith's  Pilot  Model 


Figure  4  shows  a  typical  •  vs  u  curve  (the 

“z°z 

computation  of  *a  a  is  discussed  later).  The  center 

frequency  UR)  and  the  "width"  (a')  are  the  important 

parameters  for  PIO  analysis.  A'  can  be  computed  as 

A '  =  02  /2H  [12] 

a  7 

zp 

where 

o2  =  1/(2  d«a  [13] 

ci^p  “  azaz 
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Figure  4.  Normal  Acceleration  Power  Spectral  Density 

The  integration  can  easily  be  carried  out  numerically 

with  reasonable  values  for  the  limits,  since  ♦  usually 

azaz 

tends  toward  zero  after  a  finite  value  of  frequency  'if  it 

doesn't,  *  is  certainly  not  narrowband).  From  the  above 
azaz 

computations,  an  "index  of  subjective  predictability"  (v)  can 
be  defined,  where 

v  =  A  '  /<i>p  [14] 

Reference  10  suggests  that  an  input  stimulus  to  a 
pilot-vehicle  system  may  be  "subjectively  predictable"  when 
v<0.3;  i.e.,  the  pilot  may  at  some  point  switch  to  tracking 
normal  acceleration  rather  than  pitch.  This  is  the  basis 
behind  the  occurrence  of  Type  I  PIO. 
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Staith  postulates  that  Type  I  PIO  begins  with  highly 

resonant  closed  loop  pitch  attitude  dynamics.  If  a 

pilot-vehicle  system  is  nonresonant  for  a  reasonable  variety 

of  turbulence  and  pilot  model  parameters,  then  Type  I  PIO  is 

unlikely  for  that  configuration.  During  normal  flight,  the 

pilot  model  shown  in  Figure  3  is  used  with  the  "switch"  set 

to  prevent  a  feedback  control .  The  servo  pilot  model  is 
Zp 

then  combined  with  the  aircraft/control  system  dynamics  to 

determine  if  the  pilot-vehicle  system  is  nonresonant.  If 

*  (the  acceleration  response  power  spectral  density  due 
azaz 

only  to  closed  loop  control  of  pitch  attitude)  is  determined 
to  be  sufficiently  narrowband  as  discussed  previously,  then 
Type  I  PIO  is  a  possibility. 

The  servo  model  used  in  this  study  is  the  basic  model 
from  Reference  9.  The  frequency  domain  definition  of  the 
servo  model  is  given  by 


K  e"juT(T  ju+l) 

y  -  _ P _ " _ 

P  (TNjw+l)(TIju+l) 

where 


Kp  =  pilot  gain 
t  =  reaction  time  delay 
(TLj»+l) 

-  =  equalization  characteristics 

(T  j w+1 ) 

1 

-  =  neuromuscular  system  characteristic 

(Tjgj  U+1 ) 


[15] 
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Tn  is  ti.e  first  order  lag  term  included  to  describe  the 
neuranuscular  characteristics.  It  is  often  combined  with  t 
to  give  an  equivalent  Te  =  T  +  TN»  since  it  is  difficult  to 
distinguish  between  the  two.  x  varies  from  about  0.1  to  0.2 
seconds  and  from  0.1  to  0.6  seconds,  depending  upon  the 
task  and  the  physical  state  of  the  pilot.  A  nominal  value  of 
xe=0.3  seconds  has  been  successfully  used  in  many  studies  and 
will  be  used  here  without  further  comment. 

The  values  of  TL  and  Tj  are  selected  by  the  pilot 
depending  on  the  controlled  element.  Figure  5  shows  the 
compensatory  form  selected  by  the  pilot  for  aircraft  dynamics 
of  the  form  K,  K/s,  and  K/s  .  In  general,  the  pilot  chooses 
his  compensation  to  provide  an  overall  open  loop 
pilot-vehicle  system  with  K/s  type  dynamics  in  the  region  of 
the  crossover  frequency.  There  is  good  evidence  that  T^  =  0 
is  often  a  valid  model;  i.e.,  a  lead-only  model  (11).  Pilots 
in  general  do  not  like  to  generate  lag  terns  and  will  do  so 
only  if  absolutely  necessary  to  get  good  low  frequency 
response.  Lag  terms  are  usually  generated  only  for 
aircraft/ flight  control  system  combinations  which  have  a 
relatively  flat  amplitude  response  at  low  frequencies. 
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Controlled  Compensation 
Element  Chosen  Comments 


K 


K/s 

K/s2 


K  (Tt ju+l)e-Te^u  Chosen  only  if  necessary  for 

_E _  good  low  frequency  response. 

(T  ju+1)  Often,  however,  destabilizing 

1  effects  introduced  by  the  lag 

must  be  overcome  at  higher 
frequencies  by  a  lead. 


Kp(TLju+l)e“Te^ 


Generally  chosen  to  partially 
compensate  for  reaction  time 
and  neuromuscular  lags  or  in 
an  attempt  to  increase  system 
crossover  frequency. 


Figure  5.  Pilot  Compensation  for  Common  Controlled  Elements 


With  a  first  order  Pade'  approximation,  the  lead-only 
pilot  model  becomes 


(0. 5t  jw-1 ) 

(T  ju+l) - - - 

L  (0.5Teju+l) 


[16] 


The  model  can  be  parameterized  by  knowing  the 
approximate  crossover  frequency  the  pilot  chooses  for  the 
overall  pilot/system  dynamics;  i.e.,  where  the  pilot  chooses 
to  close  the  loop.  An  estimate  of  <i>c  is  given  in  Reference  5 
as 


wc  =  6.0  +  0.24m 

m  =  l/5[-|e/Fs( jl.O) |-|e/Fs(jl.5) |+|e/Fs( j4.0) I  [17] 
+  |6/F  ( j5 .0) |- j  0/F  ( j 2 . 5 ) 1+ | e/F  ( j6.0) | 

S  5  o 
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where  "m"  is  the  average  slope  of  the  magnitude  of  the  e/F 

s 

transfer  function  (in  dB/octave)  over  the  region  of  1.0  to 
6.0  radian/second. 

A  good  representative  value  of  is  0.5  seconds  (8,12). 
Two  other  values  of  TL»  1.5  and  2.5  seconds,  were  used 
initially  in  the  study  as  well,  but  changing  the  lead  to 
these  higher  values  had  no  significant  impact  on  the  results. 
Hence,  T=0.5  seconds  was  used  throughout  the  study. 

Using  the  values  just  discussed  for  the  two  parameters, 
the  pilot  model  becomes 


(0.5t  j tu-1 ) 

Y  =  K  U  )(0.5j«+l) - fi - 

P  ^  ( 0 . 5  t  ©  j  u»+l ) 


[18] 


This  model  is  plotted  with  the  e/Fg  dynamics  (Yp[e/Fg]) 

and  K  is  chosen  to  give  the  desired  crossover  frequency,  u>  . 

P  c 

As  was  previously  noted,  satisfactory  results  cannot  be 
achieved  for  some  types  of  configurations  using  the  lead-only 
pilot  model.  Configurations  which  exhibit  a  fairly  flat 
amplitude  response  at  low  frequencies  require  a  pilot  model 
which  provides  lag  compensation  (provided  the  lag 
compensation  does  not  make  the  aircraft/pilot  system  unstable 
at  the  desired  crossover  frequency).  The  NT-33  airframe 
pitch  loop  dynamics  sometimes  resulted  in  a  configuration 
with  a  fairly  flat  pitch  loop  amplitude  response  in  the 
region  of  0.5  to  2.0  radian/second.  Figure  6  depicts  two 
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Figure  6.  Example  of  Classical  Aircraft  Response  versus  "Flat"  Aircraft  Response  (Pitch  to  Stick  Force) 


LAHOS  configurations,  one  with  a  flat  region  and  one  without. 
For  configurations  with  a  flat  amplitude  response,  a  lead-lag 
servo  pilot  model  was  used,  parameterized  as  follows: 


Vj“) 


K  (T.  ju+l)e’TeJ“ 
P  L  _ 

(Tjju+1) 


[19] 


with  Tl  =  0.5  seconds,  Tj  *  1.4  seconds,  *  0.3  seconds, 
and  Kp  selected  as  in  equation  [17].  The  value  of  was 
selected  to  approximately  cancel  the  pitch  loop  zero  which 
contibuted  to  the  flat  response. 

Once  the  pilot  model  is  selected,  the  closed  loop 
dynamics  are  computed  as: 


[Y_(j-)][«e/F.(j*)][e/«e( j»>3 

e/ e  ( j ai )  =  - £ - ? -  [20] 

U  +  [Y  (j«)][6e/Fs(j»)][e/«e(ju]> 

6 c ( t )  is  an  equivalent  command  input  to  Figure  3  due  to 
vertical  w-gusts: 


ec(  j  oi )  =  [-e/wg(  ju)]wg(  ju) 


[21] 


The  gust  response  transfer  function  0/wg  was  derived  in 


equation  [10]. 
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In  order  to  determine  the  gust  response,  a  turbulence 
model  must  be  selected.  For  this  study,  the  Dryden 
turbulence  power  spectral  density  fran  Reference  3  was 
chosen : 


w  w 

g  g 


(*>)  = 


o^L  /* 
w  w' 


i+3(y./vo)2 

[l+(Lwu,/Vo)2]2 


1 22] 


where 

VQ  =  reference  velocity  of  aircraft 

ow  =  "strength"  of  turbulence  in  feet/second 

L  =  "scale"  of  turbulence  in  feet 
w 

From  Reference  13,  ow  =  5  feet/second  is  representative 
of  moderate  turbulence,  and  this  value  was  used  throughout 
the  study.  A  value  of  Lw  =  50  feet  is  sufficient  to 
represent  the  landing  phase  of  flight. 

Now  the  power  spectral  density  of  the  pilot's  normal 
acceleration  due  to  pitch  attitude  tracking  can  be  derived. 
The  closed  loop  pitch  power  spectral  density  is 


=  I  0/w_(j<*>)  1 2*  w  U)  [23] 

Me  9  wgwg 

The  corresponding  normal  acceleration  due  to  closed  loop 
pitch  attitude  tracking  is  then 


[N^p(j.)/N8  (j. )][•/«  <»)) 

6  e  6  e  C 


[24] 
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and  finally,  the  power  spectral  density  of  the  normal 
acceleration  due  to  closed  loop  control  of  pitch  attitude  can 
be  obtained: 


,azaz(*,)  =  lazp/#c<3“')  I  V.CU) 


[25] 


where  the  pilot's  normal  acceleration  dynamics  were  derived 
previously. 

The  index  of  subjective  predictability  (v)  can  then  be 

calculated.  If  v>0.3  or  if  ♦  has  no  prominent  center 

azaz 

frequency  for  all  choices  of  pilot  and  turbulence  models, 
then  Type  I  PIO  is  unlikely.  If  more  than  one  center 
frequency  exists,  then  all  of  those  frequencies  for  which 
v^0.3  must  be  used  to  evaluate  the  phase  and  magnitude 
criteria  presented  below.  If  the  center  frequency  meets  the 
phase  and  magnitude  criteria,  then  it  should  be  considered  a 
potential  PIO  frequency. 

At  this  point,  Sftiith  suggests  that  a  simpler  criterion 

for  subjective  predictability  may  be  based  on  the  dominant 

mode  damping  ratio  cCL  for  closed  loop  pitch  attitude 

control.  He  suggests  that  a_  (t)  is  subjectively  predictable 

Zp 

when  c_  <0.2  for  the  closed  loop  pitch  to  stick  force 
CL/ — 

dynamics.  The  corresponding  center  frequency  should  be  set 

to  the  dominant  mode's  undamped  natural  frequency.  Although 

the  more  detailed  method  of  determining  subjective 

predictability  was  used  throughout  this  study,  the  ;_T < 0.2 

CL/— 


criterion  was  evaluated  as  well. 
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If  the  closed  loop  control  of  pitch  attitude  produces  a 

«  which  is  subjectively  predictable,  then  Snith 
azaz 

postulates  that  the  pilot  will  attempt  to  track  a^  at  sane 

point;  i.e.  the  pilot  behaves  as  in  Figure  3  with  the  switch 

set  to  prevent  pitch  control.  He  further  postulates  that 

when  the  pilot  begins  to  track  a  ,  his  crossover  frequency 

ZP 

will  equal  the  resonant  frequency  of  the  a  (jw)  response  due 

ZP 

to  pitch  attitude  control,  and  K  will  be  selected  to 

d 

establish  this  crossover  frequency.  The  pilot  model  then 
becomes 


F  /a  (jw)  =  K  e-Ta^u  [26] 

s  Zp  a 

with  t  =0.25  seconds.  This  value  of  t  was  chosen  by  Smith 
a  a 

because  of  its  consistency  with  actual  PIO  experiences. 

If  no  important  nonlinearities  exist  (as  assumed  in  this 
study)  then  a  necessary  condition  for  Type  I  PIO  is  that  the 
phase  margin  of  the  aZp/Fs  system  at  the  resonant  frequency 
resulting  fran  e/F  closure  must  be  negative;  i.e., 

♦  =  180°+$(jw_)  <  0, 

m  k 

where  $(jw  )  is  the  phase  angle  of  a  /F  (jw  ). 

K  Z  p  S  K 

If  the  phase  margin  is  positive,  then  Type  I  PIO  cannot 
occur,  because  a  switch  to  tracking  normal  acceleration  will 
result  in  a  stable  system.  To  determine  the  phase  margin, 
♦(jw)  is  plotted  for  the  open  loop  a  /F  dynamics  including 

Z  p  s 

a  pilot  phase  resulting  from  a  0.25  second  delay,  as  in 


equation  [25].  ♦(jw)  will  be  the  sum  of  phase  angles  due  to 

the  pilot,  the  feel  system  dynamics,  the  control  system 
dynamics,  and  the  airframe  dynamics. 

Smith’s  final  necessary  condition  for  Type  I  PIO  is  that 

the  amplitude  of  a  oscillations  due  to  pitch  attitude 

ZP 

control  must  be  greater  than  seme  "critical"  value  (a  )  to 

Z p  CR 

cause  the  pilot  to  make  the  "switch"  to  acceleration 
tracking.  It  is  postulated  that  (a  )_  =  0.012 

Zp  tK 

g/degree/second,  i.e. 


lazp/®(  j“R>  I 


>  0.012  g/degree/second 


is  necessary  for  the  pilot  to  sense  a  subjectively 
predictable  normal  acceleration  power  spectrum  due  to  pitch 
attitude  tracking.  This  condition  will  be  called  "Smith's 
magnitude  criterion." 

In  summary,  for  Type  I  PIO  with  no  significant 
nonlinearities : 

1.  Select  an  appropriate  pilot  model  for  the 
aircraft/ flight  control  system  dynamics  based  on  the  form  of 
the  dynamics  and  the  crossover  frequency  computation 
described  in  Equation  [17]. 

2.  Close  the  pitch  attitude  loop  (e/F  ). 

s 

3.  Compute  a  ,  (w)  using  a  representative  Dryden  model 

a  z«  z 

of  vertical  turbulence. 
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4.  Estimate  «R  from  #a  a  •  If  <*>R  exists  for  any 
pilot/vehicle  combination  then  estimate  the  subjective 
predictability  index  v  If  v>0.3  then  go  to  step  4a; 
otherwise,  go  to  step  5. 

4a.  Estimate  the  resonant  mode  damping  ratio  tR.  If 

C  >0.2  then  Type  I  PIO  is  unlikely.  If  C  <0.2  go  to  step  5. 

R  R 

5.  Plot  the  total  open  loop  system  phase  angle  Bode 
♦(jw)  for  the  &  /F  loop  dynamics,  ♦(jw)  will  be  the  sum  of 

Z  p  S 

phase  angles  due  to  the  pilot,  the  feel  system  dynamics,  the 
control  system  dynamics,  and  the  airframe  dynamics.  The 
pilot  phase  should  be  assumed  to  result  entirely  frcm  a  0.25 
second  delay.  If  the  phase  margin  180°+*(ju  )<0  then  go  to 
step  6. 

6.  If  la,,  / e  ( jwD)  I  <0.012  g/degree/second  then  conclude 

Z  p  A 

that  Type  I  PIO  is  unlikely.  If  this  ratio  is  >  0.012,  then 
conclude  that  Type  I  PIO  is  a  possibility. 

Type  II  PIO.  Type  II  PIO  is  initiated  when  an  abrupt 
control  or  disturbance  of  sufficient  amplitude  excites  the 
stick-free  dynamic  modes  of  the  aircraft.  This  type  of 
control  might  be  due  to  open  loop  high  g  maneuvering,  system 
transients  from  SAS/CAS  shutdown  or  start-up,  and  so  on. 
Although  by  definition  Type  II  PIO  is  not  likely  to  occur 
during  the  landing  phase  of  flight,  it  is  presented  here  for 
completeness  of  smith's  theory. 

Type  II  PIO  is  analyzed  similarly  to  Type  I  PIO,  except 
that  only  open  loop  dynamics  are  considered.  The  procedure 
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is  essentially  the  same,  but  there  is  no  loop  closure  to 

determine  ♦  .  Instead,  ♦  is  computed  using  a 

«zaz  «zaz 

normalized,  broadband  noise  representation  for  Fg(ju)  to 
simulate  the  required  "abrupt"  character,  i.e. 


* 

a  za 


( w) 

z 


|a  /F  (s)| 
ZP  8 


2 

s=j  u 


1 


The  procedure  then  follows  the  same  steps  as  Type  I  PIO. 
smith  also  proposes  that  a  simpler  criterion  for  Type  II 
PIO  may  be  to  determine  whether  any  stick-free  dynamic  mode 
exists  with  damping  ratio  ;R_<0.2  which  significantly 


contributes  to  ♦  _  (w).  The  criterion  for  "significant 

azaz 

contribution"  is  that  the  modal  frequency  is  less  than  about 
10  radian/ second  (the  pilot's  bandwidth  of  control  limit). 
Again,  if  such  a  mode  exists,  the  response  is  subjectively 
predictable,  <uR  is  then  set  to  the  modal  frequency,  and  the 
additional  criteria  are  examined  at  wR. 

In  summary,  for  Type  II  PIO, 

1.  Compute  the  power  spectral  density  (PSD)  of  a  for 

ZP 

the  stick-free  airplane  dynamics;  assume  that  the  airplane 
is  excited  by  a  wideband  noise  with  PSD=1  (to  simulate  abrupt 
inputs).  That  is,  assume 


azaz 


U)  =  |a  /F  (s)  |f  .«  x  1 
zD  s  s=j 


2.  Continue  with  the  analysis  as  described  above  for 
Type  I  PIO,  starting  at  step  4  and  replacing  "Type  I"  with 
"Type  II." 
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3.  Simplified  alternative  procedure:  if  the  damping 

ratio  of  the  dominant,  resonant  mode  of  a  /F  is  <0.2,  then 

zp'  s  - 

Type  II  PIO  is  possible.  Call  this  damping  ratio  ;  and 
continue  the  analysis  as  described  above  for  Type  I  PIO, 
starting  with  Step  5  and  replacing1  "Type  I"  with  "Type  II." 

If  the  damping  ratio  is  >0.2,  then  conclude  that  Type  II  PIO 
is  unlikely.  For  conservatism,  one  could  define  wR  as  the 
dominant  mode's  damped  frequency  and  proceed  to  step  5  of  the 
Type  I  PIO  analysis  above. 

Type  III  PIO.  Type  III  PIO  of  Reference  5  is  postulated 

to  appear  during  the  single  loop  tracking  of  pitch  attitude; 

i.e.,  the  pilot's  normal  acceleration  dynamics  are 

irrelevant.  This  mode  will  probably  only  be  seen  when 

control  system  or  equivalent  time  delays  induce  significant 

phase  lag  within  the  bandwidth  of  the  pilot's  control  (about 

10  radian/ second ) .  The  time  delay  must  be  sufficient  to  make 

the  e/F  loop  unstable  at  w  (the  crossover  frequency 
s  c 

predicted  by  Equation  [17]);  i.e.,  the  necessary  condition 

for  attitude  only  PIO  is 

N/F  (ju  )  <  -180° 
s  c 

To  get  PIO,  something  has  to  excite  the  pilot  to  adapt 
the  form  of  a  pure  gain.  Smith  assumes  that  this  follows  the 
development  of  substantial  aircraft  resonance  in  attitude 
response,  as  for  Type  I  PIO,  except  now  the  pilot's  time 
delay,  shown  in  Figure  3,  is  assumed  to  be  zero  (this  seems 
to  be  for  lack  of  a  better  value).  The  pilot  also  either 
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feels  no  normal  acceleration  (as  in  a  fixed-base  simulator) 
or  does  not  attempt  to  track  the  normal  acceleration  he 
feels. 

The  technique  for  Type  III  PIO  analysis  is  to  use  a  pure 
gain  pilot  model  and  look  for  locus  crossings  of  the  ju-axis 
on  a  root  locus  plot.  When  the  frequency  at  axis  crossing  is 
less  than  u>c,  then  the  crossing  conditions  (the  pilot  gain  at 
the  crossing  and  the  corresponding  resonant  frequency) 
represent  a  potential  PIO  state;  i.e.,  the  pilot  may  adopt 
the  pure  gain  model  and  cause  PIO. 

In  summary,  for  Type  III  PIO: 

1.  Plot  the  root  locus  of  e/F  . 

s 

2.  If  there  are  no  imaginary  axis  crossings  of  the  root 
locus,  then  Type  III  PIO  is  not  a  possibility. 

3.  If  the  root  locus  crosses  the  imaginary  axis, 
determine  the  frequency  of  crossing  (wD).  If  iuD<u  ,  where  u 
is  determined  from  Equation  [17],  then  Type  III  PIO  is  a 
possibility.  The  gain  of  the  root  locus  at  the  axis  crossing 
determines  the  gain  the  pilot  must  generate  to  cause  PIO. 

Bandwidth  Method 

The  bandwidth  method  of  Reference  4  was  proposed  as  a 

handling  qualities  requirement  for  MIL-F-8785C.  The 

bandwidth  method  is  particularly  attractive  because  it 

assumes  a  "gain-only"  pilot  model  and  involves  only  the  use 

of  open  loop  pitch  to  stick  force  (e/F  )  Bode  plots. 

s 
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"Bandwidth"  can  be  loosely  defined  as  the  maximum 
frequency  at  which  closed  loop  compensatory  tracking  can  take 
place  without  threatening  the  stability  of  the  aircraft;  i.e. 
the  maximum  open  loop  crossover  frequency.  Hence,  a  large 
value  of  bandwidth  is  generally  desirable  to  achieve  superior 
tracking  performance. 

The  reason  for  including  bandwidth  in  this  study  was  to 
determine  if  the  simple  bandwidth  method  alone  could  be  used 
to  separate  PIO  prone  aircraft  from  those  which  are  not.  The 
bandwidth  criterion  sets  up  boundaries  for  Level  1,  2,  and  3 
handling  qualities  for  both  Category  A  and  C  requirements, 
based  on  maximum  crossover  frequency  and  system  phase  delay. 
The  approach  taken  in  this  study  was  to  see  how  well  the 
boundaries  and  the  phase  delay  parameter  correlated  with  PIO 
tendencies  noted  during  the  LAHOS  study. 

The  following  discussion  of  the  bandwidth  theory  is 
taken  from  reference  4.  Crossover  frequency,  directly 
determined  by  pilot  gain,  is  a  rough  measure  of  the  rapidity 
of  a  closed  loop  response.  Physically,  the  pilot  will 
increase  his  gain  (and  hence,  crossover  frequency)  to  track 
more  rapidly  moving  targets  with  acceptable  error.  However, 
the  pilot  cannot  indefinitely  increase  crossover  frequency  by 
increasing  gain,  because  he  will  eventually  lose  closed  loop 
stability  (when  the  phase  margin  of  the  open  loop  system 
becomes  negative)  .  The  pilot  would  like  to  choose  a  value 
which  allows  him  to  double  his  gain  and  provide  plently  of 
phase  margin.  A  reasonable  crossover  frequency  would  then  be 


33 


one  which  provides  at  least  6  dB  of  gain  margin  and  45°  of 
phase  margin. 


The  above  crossover  frequency  is  the  bandwidth  frequency 
UBW)  and  is  shown  in  Figure  7.  uQW  is  defined  to  be  the 
smaller  of  two  values,  »BWphase  and  *BWphase.  These  two 
frequencies  are  determined  as  shown  in  Figure  7. 


Frequency  Response  of  Pitch  to 
Stick  force  Dynamics 


Frequency  (red/sec) 

Tp  =  -(;{,2(i)1  80  +  180o)/(57.3x2(d180) 


Figure  7.  Definition  of  Bandwidth  and  Phase  Delay  Parameters 


Handling  qualities  and  pilot  ratings  are  not  dependent 

on  bandwidth  alone;  the  shape  of  the  phase  curve  at 

frequencies  above  u>fiW  becomes  important  as  well.  If  the 

phase  curve  drops  off  rapidly  at  frequencies  above  u  ,  the 

BW 

aircraft  will  generally  receive  poor  pilot  ratings,  since  an 
abrupt  loss  in  stability  margin  is  produced  when  the  pilot 
attempts  to  increase  the  crossover  frequency.  One  measure  of 
rapid  phase  rolloff  is  equivalent  system  time  delay. 

However,  equivalent  system  time  delay,  unlike  uaw»  is  not 
easily  measured.  Phase  delay,  x  ,  a  parameter  which  is 
easily  measured,  is  defined  in  Figure  7.  Usually,  is 
numerically  similar  to  equivalent  time  delay. 

The  bandwidth  criterion  suggests  that  systems  with  high 
attainable  crossover  frequencies  and  without  rapid  phase 
rolloffs  should  have  good  handling  qualities.  Figure  8,  from 
Reference  4,  shows  flying  qualities  boundaries  based  on 
bandwidth  frequency  and  x  .  The  boundaries  of  Figure  8  are 
referred  to  as  the  "bandwidth  criterion." 
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(sec) 


Category  A  Flight  Phases  Category  C  Flight  Phases 


Figure  8.  Hoh's  Bandwidth  Criterion  for  Handling  Qualities 

The  bandwidth  method  can  then  be  summarized  as.  follows: 

1.  Determine  u>1Q0  from  the  Bode  plot  of  e/Fg. 

2.  Find  “BWphase=frequency  where  the  phase  margin  is 

45  . 

3.  Find  | 9/Fg |  at  wiao  and  add  6dB.  Call  this  value  Z. 

4.  Find  uBWgain=f requency  where  Z  occurs. 

5‘  UBW  =  min((UBWphase'  “BWgain^ 

6‘  TP  =  ^ ( 2u180  )+180°]/[ 57 . 3x2u>180]  . 

7.  Find  Tp  and  wBW  on  the  vs  wBW  plot  to  determine 
the  predicted  level  of  handling  qualities. 


YF-17  Example 


An  example  will  now  be  presented  which  applies  the  two 
theories  to  the  YF-17  as  simulated  using  the  NT-33  during  the 
LAHOS  study.  The  original  YF-17  (as  simulated  on  the  NT-33) 
sustained  a  severe  PIO  during  the  landing  flare;  a  modified 
version  was  PIO  free. 

The  transfer  functions  for  the  YF-17  ares 
8  0.3369(0.0853) (0.6870) 

-  SS  ••  '  - 

«e  [0.15,0.16][0.65,1.94] 

where 

(a)  =  (s+a) 

2  2 

[  c,  u  ]  =  (s  +2c<*»  S+w  ) 
n  n  n 

ar  -1.066(0)(0. 0258)[-0.02,6. 80] 

..  -Z..E-  =  - 

«e  [0. 15,0. 16][0. 65, 1.94] 

6  -0.0085(0) (-0. 0051) (-0.06528) 

wg  [0. 15,0. 16][0. 65, 1.94] 

5625.0 

Actuator  Dynamics:  - 

[0.07,75.0] 


Feel  System: 


84.5 

[0.6, 26.0] 


16.37(2.0) (2.3) 

Unmodified  - 

Control  System:  -  (0.9)(5*0){0.7,4.03 

4.26(2.0) (2. 3)  (16. 7) 

Modified  - 

Control  System:  ( 0 . 9 ) ( 5 . 0 ) ( 10. 0) 

! 
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Smith's  Theory.  Tables  1  and  2  summarize  the  results  of 
applying  the  Type  I  PIO  theory  to  the  YF-17  data.  Table  1 
shows  the  results  using  the  v_<0.3  criterion  and  Table  2  shows 
the  results  using  the  t<0.2  criterion. 


Table  1 


Results  of  Type  I  PIO,  YF-17, 

Using  Index  of  Subjective  Predictablity 


Aircraft  wR(rad/sec) 

V 

Margin(deg) 

(g/deg/sec) 

Original  0.17 

125.88 

. 

YF-17  0.85 

0.13 

104.76 

— 

2.69 

0.03 

-40.61 

0.0241 

9.24 

1283.18 

- 

- 

Modified  0.17 

23.62 

_ 

YF-17  4.22 

0.06 

-38.22 

0.0114 

^Phase  Margin  =  phase  margin  of 

a  /F  evaluated 

Zr,  8 

at  uR  using 

pure  gain  plus  time  delay  (t  =0. 

a 

r 

25  seconds)  for  pilot  model 

Results 

Using  Dominant 

Table  2 

of  Type  I  PIO,  YF-17, 

Closed  Loop  Resonant  Frequency 

Aircraft  uR( rad/sec ) 

Ccl 

Phase 

Margin(deg) 

(g/deg/sec) 

Oriqinal  0.84 

-0.L2  _ _ 

105.63 

2.73 

-0.08 

-43.05 

0.0236 

Modified  4.60 

0.14 

-49.10 

0.0092 
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The  original  YF-17  has  two  subjectively  predictable 

resonant  frequencies,  one  at  about  0.85  radian/ second,  and 

another  at  about  2.7  radian/second.  The  frequencies 

predicted  by  both  the  vjc0.3  criterion  and  the  ccl_50.2 

criterion  are  very  close,  a  /F  has  a  positive  phase  margin 

Zp  s 

at  u  =0.85  radian/ second,  but  is  negative  at  u  =2.7 
R  R 

radian/ second.  The  magnitude  criterion  at  this  frequency  is 

about  twice  the  value  needed  for  PIO  when  evaluated  at  w  =2.7 

R 

radian/ second.  Hence,  the  original  YF-17  should  have  a  high 

probability  of  being  prone  to  Type  I  PIO. 

The  modified  YF-17  has  a  subjectively  predictable 

resonant  frequency  at  uir=4.22  radian/second  using  the  v<0.3 

criterion.  The  Ccl_f0.2  criterion  predicts  a  frequency  of  4.6 

radian/second .  These  values  are  within  10  percent  of  each 

other.  Using  either  frequency,  the  phase  margin  criterion  is 

negative,  but  the  magnitude  criterion  is  not  met.  Therefore, 

Type  I  PIO  is  unlikely  for  the  modified  YF-17. 

Table  3  summarizes  the  results  of  applying  the  Type  II 

PIO  theory  to  the  YF-17  data.  This  table  summarizes  only  the 

v_<0.3  criterion,  since  no  frequencies  were  predicted  using 

the  c  <0.2  criterion, 
cl— 
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Table  3 


Results  of  Type  II  PIO,  YF-17, 

Using  Index  of  Subjective  Predictability 


Phase 


|az  /•(■„)  I 


Aircraft 

wR(rad/sec) 

V 

Margin (deg) 

(g/deg/sec) 

Original 

0.17 

0.09 

225.91 

- 

YF-17 

9.91 

371.34 

- 

- 

Modified 

0.17 

0.10 

228.92 

Phase  Margin  =  phase  margin  of  a  /F  evaluated  at  wD  using 

Zp  S  K 

pure  gain  plus  time  delay  (t  *0.25  seconds)  for  pilot  model 


Although  each  configuration  has  a  subjectively 
predictable  resonant  frequency  at  0.17  radian/ second,  both 
have  large  positive  phase  margins.  Hence,  Type  II  PIO  is 
unlikely  for  either  configuration. 

The  root  loci  of  both  the  original  and  modified  YF-17 
are  shown  in  Figure  9.  The  jw-axis  crossing  is  annoted  on 
each  diagram.  The  predicted  crossover  frequency  from 
equation  [17]  and  the  jw-axis  crossing  are  compared  in  Table 
4. 
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-12.0  -10.0  -8.0  -6.0  -4.0  -2.0  0.0  2J 

REAL  AXIS 
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REAL  AXIS 


Figure  9.  Root  Loci  for  Original  and  Modified  YF-17 


Table  4 


Predicted  Crossover  Frequency  and  jw-axis  Crossing,  YF-17 


Aircraft 

Predicted  u  * 
(rad/sec)c 

ju-axis  crossing 
(rad/sec) 

Original 

3.18 

2.34 

YF-17 

Modified 

YF-17 

3.84 

5.86 

^Predicted  wc  =  crossover  frequency  predicted  by  Smith's 
formula,  based  on  average  slope  of  the  pitch  to  stick 
force  magnitude  curve 

Type  III  PIO  would  be  unlikely  for  the  modified  YF-17. 
For  the  modified  YF-17,  the  frequency  at  the  jw-axis  crossing 
is  less  than  the  predicted  crossover  frequency.  Hence,  the 
unmodified  YF-17  should  also  be  suspected  of  being  prone  to 
Type  III  PIO  (caused  by  excessive  time  delay  in  the  system). 
How  to  determine  which  type  of  PIO  was  actually  experienced 
during  the  flight  test  is  unclear,  but  it  may  be  related  to 
the  actual  PIO  frequency.  Type  I  PIO  theory  predicts  a  PIO 
frequency  of  about  2.7  radian/ second;  Type  III  PIO  theory 
predicts  a  PIO  frequency  of  about  2.3  radian/second.  The 
actual  PIO  value  was  about  3.0  radian/second;  hence,  the 
original  YF-17  configuration  probably  experienced  a  Type  I 
PIO. 
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Bandwidth  Method.  Bode  plots  for  both  the  original  and 
modified  YF-17  are  shown  in  Figures  10  and  11.  The  critical 
values  needed  for  the  bandwidth  method  are  shown  on  the 
diagrams.  Figure  12  depicts  the  bandwidth  criterion  for  both 
configurations.  The  modified  YF-17  (with  no  PIO  tendency)  is 
very  close  to  the  Level  1  boundary  and  is  predicted  to  have 
good  handling  qualities.  The  original  YF-17  (which 
experienced  a  PIO)  is  predicted  to  have  Level  3  handling 
qualities.  This  single  data  point  suggests  the  bandwidth 
criterion  boundaries  might  be  useful  as  PIO  predictors. 
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Figure  10.  Critical  Values  used  for  Bandwidth  Method,  Original  YF-17 
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Figure  12.  Bandwidth  Criterion  for 

Original  and  Modified  YF-17 
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III.  Analytical  Results  and  Analysis 


To  perform  a  preliminary  analysis,  the  techniques  were 
applied  to  the  Calspan  LAHOS  data  base  (6).  This  data  base 
was  the  result  of  a  1978  Calspan  study  using  the  variable 
stability  NT-33A  to  determine  the  effects  of  higher  order 
control  system  dynamics  on  fighter  approach  and  landing 
handling  qualities.  This  was  an  ideal  data  base  to  use 
because  many  PlOs  were  experienced  during  the  flight  tests. 
Appendix  A  contains  a  complete  description  of  the  LAHOS 
configurations . 

Before  beginning  an  analysis  of  the  data,  it  is 
necessary  first  to  define  what  is  meant  by  a  "PIO  tendency" 
and  develop  a  method  to  determine  if  a  PIO  tendency  has  been 
experienced  in  flight.  Figure  13  contains  both  the  PIO 
rating  scale  and  a  flowchart  which  is  commonly  used  in-flight 
by  pilots  to  determine  PIO  rating.  From  the  definitions  on 
the  scale,  there  is  no  doubt  that  a  rating  of  4,  5,  or  6  can 
be  classified  as  a  PIO  tendency.  Likewise,  a  rating  of  1 
indicates  no  PIO  tendency.  The  question  then  arises  as  to 
whether  or  not  the  "undesirable  motions"  necessary  for  a  PIO 
rating  of  2  or  3  ought  to  be  considered  a  PIO  tendency.  A 
full-blown  PIO,  by  definition,  implies  the  existence  of  a 
sustained  oscillation;  however,  undesirable  motions  such  as  a 
bobble  or  a  tendency  to  overcontrol  usually  precede  a  PIO. 
What  prevents  these  undesirable  motions  from  becoming  PIOs  is 
that  the  pilot  is  able  to  provide  enough  compensation  to  keep 
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the  motion  down  to  a  nuisance  level.  Thus,  for  this  study,  a 
PIO  tendency  was  defined  to  be  a  PIO  rating  of  2  or  greater. 

Simply  defining  a  PIO  tendency  is  not  enough  to 
determine  if  a  PIO  tendency  actually  existed,  however.  If  a 
configuration  is  flown  several  times,  it  is  possible  that 
ratings  may  be  averaged,  with  an  average  PIO  rating  2  or 
greater  indicating  a  PIO  frequency. 

There  are  two  problems  associated  with  averaging  PIO 
ratings.  The  first  problem  is  that  the  PIO  rating  scale  is 
not  necessarily  linear  in  the  sense  that  the  severity  of  the 
PIO  may  not  increase  linearly  with  increasing  PIO  rating. 
Since  there  is  no  present  means  of  quantifying  PIO  severity, 
there  is  no  way  to  know  if  the  relationship  is  linear  or  not. 
For  this  reason,  although  the  ratings  used  in  this  study  were 
averaged,  the  spread  of  actual  ratings  for  each  configuration 
is  included  for  completeness. 

The  second  problem  arises  with  configurations  which  are 
flown  only  once  or  twice,  as  many  of  the  LAHOS  configurations 
were.  Particularly  poor  aircraft  (aircraft  which  receive  PIO 
ratings  of  3  or  greater)  tend  to  provide  repeatable  results; 

1. e.,  virtually  everyone  who  flies  the  configuration  will 
give  it  the  same  rating.  Therefore,  even  if  a  particularly 
poor  configuration  is  flown  only  once,  there  can  be  a  fairly 
high  level  of  confidence  that  the  PIO  rating  is  correct.  The 
problem  exists  with  the  ratings  of  1  and  2.  In  general,  one 
data  point  is  insufficient  to  distinguish  between  a  1  and  a 

2.  The  LAHOS  data  had  several  configurations  which  received 
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ratings  of  1  or  2  or  both  which  were  flown  only  once  or 
twice.  Table  5  summarizes  the  number  of  flights  and  the  PIO 
ratings  for  the  LAHOS  configurations.  Each  configuration  is 
classified  as  PIO  prone,  not  PIO  prone,  or  unsure.  However, 
despite  the  uncertainty  in  sane  of  the  data,  all  the  LAHOS 
configurations  were  included  in  the  preliminary  analysis. 

It  must  also  be  pointed  out  that  the  argument  presented 
above  is  really  of  academic  interest  only,  since  an  aircraft 
which  receives  a  PIO  rating  of  2  is  unlikely  to  be  considered 
in  need  of  modification.  It  is  more  important  to  be  able  to 
determine  if  a  theory  can  separate  a  PIO  rating  of  1  or  2 
from  a  3,  since  a  PIO  rating  of  3  indicates  a  problem 
probably  requiring  correction. 


Table  5 


Summary  of  PIO  Ratings,  LAHOS  Configurations 


Configu¬ 

ration 

Number 

Plights 

Actual 

PIO  Rating^ 

v.  »~age 

' 

Consider 
PIO  Prone 

1-A 

1 

1 

UNSURE 

1-B 

1-C 

1 

2 

2 

1/1 

1 

UNSURE 

NO 

1-1 

2 

2/1 

1.5 

UNSURE 

1-2 

1 

2 

2 

UNSURE 

1-3 

4 

4/ 4/2/ 3 

3.25 

YES 

1-4 

1 

4 

4 

YES 

1-6 

2 

2/2 

2 

YES 

1-8 

1 

3 

3 

YES 

1-11 

1 

3.5 

3.5 

YES 

2-A 

2 

2/2.5 

2.25 

YES 

2-C 

4 

2/ 1/1/1 

1.25 

NO 

2-1 

3 

1/1/1 

1 

NO 

2-2 

2 

2/1 

1.5 

UNSURE 

2-3 

1 

3 

3 

YES 

2-4 

3 

3/2/1 

2 

YES 

2-6 

1 

2.5 

2.5 

YES 

2-7 

2 

3/3 

3 

YES 

2-9 

1 

3 

3 

YES 

2-10 

1 

4 

4 

YES 

2-11 

1 

3 

3 

YES 

3-C 

2 

1/1.5 

1.25 

NO 

3-1 

3 

2/3/2 

2.33 

YES 

3-2 

2 

3/3 

3 

YES 

3-3 

2 

4/3.5 

3.75 

YES 

3-6 

2 

3/3 

3 

YES 

3-7 

1 

4 

4 

YES 

4-C 

2 

1.5/2 

1.75 

UNSURE 

4-1 

1 

1 

1 

UNSURE 

4-3 

3 

2/3/3 

2.67 

YES 

4-4 

3 

3/3/2 

2.67 

YES 

4-6 

1 

2 

2 

UNSURE 

4-7 

1 

1 

1 

UNSURE 

4-10 

1 

4 

4 

YES 

4-11 

1 

4 

4 

YES 

5-1 

2 

3/3 

3 

YES 

5-3 

5 

3/3/2. 5/1/1 

2.83 

YES 

5-4 

1 

2.5 

2.5 

YES 

5-5 

1 

3 

3 

YES 

5-6 

1 

3 

3 

YES 

5-7 

1 

3 

3 

YES 

5-11 

1 

3.5 

3.5 

YES 

6-1 

1 

4 

4 

YES 

6-2 

1 

1 

1 

UNSURE 
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Stoith's  Theory 


To  use  Stoith's  theory,  all  three  types  of  PIO  must  be 
examined  before  concluding  whether  or  not  a  PIO  tendency  is 
predicted.  Type  II  PIO  was  not  predicted  for  any  of  the 
LAHOS  configurations,  so  it  will  not  be  included  in  the 
discussion  below. 

Tables  6  and  7  summarize  Stoith's  Type  I  PIO  theory  as 
applied  to  the  LAHOS  data.  The  lead-only  pilot  model  was 
applied  to  all  configurations,  except  as  noted.  Table  6 
summarizes  the  results  using  the  v£0*3  criterion  for 
predicting  subjectively  predictable  resonant  frequency,  and 
Table  7  summarizes  the  results  using  the  ;^0.2  criterion. 

The  two  criterion  will  hereafter  be  referred  to  as  the  "v 
criterion"  and  the  "t  criterion,"  respectively. 

The  two  criteria  are  in  very  close  agreement,  for  the 
most  part,  concerning  subjectively  predictable  resonant 
frequency  and  PIO  tendency  predictions.  However,  there  are  a 
few  notable  exceptions. 

The  v  criterion  predicts  resonant  frequencies  for 
several  configurations  which  the  C  criterion  does  not.  These 
configurations  are  summarized  in  Table  8.  All  the  extra 
frequencies,  except  one,  have  values  of  v  greater  than  0.09. 
Therefore,  vj<o.09  may  be  a  better  cutoff  for  subjective 
predictability.  The  only  configuration  which  would  have  any 
difficulty  with  this  value  of  v  would  be  4-C.  This 
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Smith's  Type  I  PIO  Predictions  Using  v  Criterion,  LAHOS 
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configuration  has  a  subjectively  predictable  resonant 
frequency  at  7.93  radian/second  with  v=0.1776.  The  C 
criterion  predicts  a  frequency  of  7.29  radian/second . 


Table  8 

Discrepancies  Between  v  and  c  Criteria 
Using  Lead-Only  Pilot  Model 


Configu¬ 

ration 

(Or 

(rad/sec) 

V 

1-A 

1.30 

0.1455 

1-B 

1.36 

0.1493 

1-C 

1.47 

0.1651 

3-1 

2.23 

0.0962 

3-2 

2.31 

0.1037 

3-6 

2.29 

0.0943 

3-7 

2.31 

0.0749 

4-C 

2.19 

0.2432 

4-1 

3.06 

0.1380 

The  other  notable  exception  is  the  frequency  difference 
predicted  by  the  two  criteria  when  the  lead/lag  pilot  model 
is  used.  In  addition,  the  5  criterion  predicted  a  resonant 
frequency  for  2-A  which  was  not  predicted  by  the  v  criterion. 
These  differences  are  summarized  in  Table  9. 
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Table  9 


Discrepancies  Between  v  and  C  Criteria 
Using  Lead/Lag  Pilot  Model 


(v<0.3)  ,  (;C1<0.2) 

Configu-  Phase1  .  Phase 


ration 

( rad/sec) 

Margin 

( rad/sec) 

Margin 

2-A 

1.08 

145.11 

1.08 

145.10 

-- 

— 

5.00 

-19.82 

2-C 

2.77 

36.92 

3.52 

6.14 

5-1 

3.71 

19.76 

4.05 

2.12 

5-6 

3.28 

25.51 

3.59 

7.92 

5-7 

3.17 

26.23 

3.45 

9.51 

5-11 

3.01 

27.80 

3.26 

12.58 

Phase  margin 

of  az  /Fs 
ZP  S 

evaluated  at 

uR,  units 

in  degrees 

The  frequency  difference  is  enough  to  make  the  phase 
margin  criterion  only  slightly  positive  for  five  of  the  six 
frequencies  predicted  by  the  C  criterion.  Using  a  slightly 
positive  phase  margin  criterion  (^15  degrees)  to  account  for 
uncertainties  in  the  lead/lag  pilot  model,  the  C  criterion 
would  then  correctly  predict  the  PIO  tendency  of  four  of  the 
above  six  configurations.  The  only  configurations  that  now 
fail  are  2-A,  which  fails  the  magnitude  criterion,  but  was 
PIO  prone  in  flight;  and  2-C,  which  is  now  predicted  to  be 
PIO  prone,  but  was  not  in  flight. 

One  other  interesting  configuration  is  4-C,  which  has  a 
very  high  subjectively  predictable  resonant  frequency  and 
very  negative  phase  margin,  but  which  fails  the  magnitude 
criterion.  There  is  reason  to  believe,  with  PIO  ratings  of 
1.5  and  2,  that  this  aircraft  may  have  had  a  PIO  tendency. 
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It  is  possible  that  the  very  negative  phase  margin  and/or 
high  frequency  may  have  caused  problems  with  Snith's  theory. 

Table  10  summarizes  Smith's  Type  III  PIO  theory.  There 
is  only  one  configuration,  5-3,  which  is  predicted  to  be  Type 
III  PIO  prone  which  was  not  also  predicted  to  be  Type  I  PIO 
prone.  Since  24  configurations  were  predicted  to  be  both 
Type  I  and  Type  III  PIO  prone,  the  obvious  question  then 
arises:  which  type  of  PIO  is  actually  seen  in  flight?  Since 
Type  III  PIO  is  dependent  only  on  excessive  time  delay  and 
not  normal  acceleration  characteristics,  a  simulator  study 
may  be  useful  in  answering  this  question. 

Table  11  summarizes  the  overall  results  of  using  Smith's 
theory  to  predict  LAHOS  PIOs.  Only  34  of  the  44  PIO 
tendencies  are  considered  fairly  certain;  of  these,  only  two, 
2-A  and  2-C  as  noted  before,  were  not  correctly  predicted. 

To  summarize  Smith's  results: 

1.  Smith's  method  works  well  as  long  as  the  lead  only 
pilot  model  is  applied.  Using  the  lead/lag  pilot  model 
produced  uncertainty  in  the  data  and  discrepancies  between 
using  the  v  and  c  criteria  to  determine  subjectively 
predictable  resonant  frequency. 

2.  There  is  little  difference  between  the  v  and  t 
criteria  for  predicting  subjectively  predictable  resonant 
frequency.  However,  the  t  criteria  seems  to  have  an  edge  on 
the  v  criterion  when  the  lead/lag  pilot  model  is  used. 
Moreover,  the  ;  criterion  is  much  simpler  to  use. 
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Table  10 


Smith's  Type  III  PIO  Predictions,  LAHOS  Data 

Configu-  “r1  ^c2  Predict 

ration  (rad/sec)  (rad/sec)  PIO? 


1-A 

11.11 

3.85 

NO 

1-B 

10.30 

3.68 

NO 

1-C 

8.44 

3.46 

NO 

1-1 

4.06 

.  3.22 

NO 

1-2 

2.41 

3.09 

YES 

1-3 

1.78 

2.74 

YES 

1-4 

1.40 

2.32 

YES 

1-6 

2.49 

3.21 

YES 

1-8 

2.06 

3.16 

YES 

1-11 

2.03 

3.12 

YES 

2-A 

12.15 

4.71 

NO 

2-C 

10.04 

4.31 

NO 

2-1 

6.48 

4.08 

NO 

2-2 

4.06 

3.95 

NO 

2-3 

3.17 

3.60 

YES 

2-4 

2.65 

3.18 

YES 

2-6 

4.12 

4.07 

NO 

2-7 

3.79 

4.06 

YES 

2-9 

2.99 

3.83 

YES 

2-10 

2.55 

3.34 

YES 

2-11 

3.44 

4.07 

YES 

3-C 

7.87 

3.86 

NO 

3-1 

3.98 

3.62 

NO 

3-2 

2.98 

3.50 

YES 

3-3 

2.59 

3.15 

YES 

3-6 

3.02 

3.62 

YES 

3-7 

2.88 

3.61 

YES 

4-C 

11.50 

4.35 

NO 

4-1 

8.21 

4.12 

NO 

4-3 

3.42 

3.64 

YES 

4-4 

2.62 

3.22 

YES 

4-6 

4.80 

4.11 

NO 

4-7 

4.30 

4.10 

NO 

4-10 

2.49 

3.38 

YES 

4-11 

3.78 

4.11 

YES 

5-1 

8.82 

5.26 

NO 

5-3 

4.66 

4.78 

YES 

5-4 

4.07 

4.37 

YES 

5-5 

3.67 

4.04 

YES 

5-6 

5.74 

5.26 

NO 

5-7 

5.28 

5.24 

NO 

5-11 

4.83 

5.26 

YES 

6-1 

2.34 

3.18 

YES 

6-2 

5.86 

3.84 

NO 

^Frequency  at  j<jj-axis  crossing 

^Crossover  frequency  predicted  using  Smith's  formula 
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Table  11 


Summary  of  Smith's  PIO  Predictions,  LAHOS  Data 


Configu¬ 

ration 

Predict 
Type  I 

Predict 
Type  III 

PIO  in 
Flight 

Actual. 
PIOR  1 

Average 

PIOR 

1-A 

NO 

NO 

UNSURE 

1 

1 

1-B 

NO 

NO 

‘ UNSURE 

2 

2 

1-C 

NO 

NO 

NO 

1/1 

1 

1-1 

YES 

NO 

UNSURE 

2/1 

1.5 

1-2 

YES 

YES 

UNSURE 

2 

2 

1-3 

YES 

YES 

YES 

4/ 4/2/ 3 

3.25 

1-4 

YES 

YES 

YES 

4 

4 

1-6 

YES* 

YES 

YES 

2/2 

2 

1-8 

YES 

YES 

YES 

3 

3 

1-11 

YES 

YES 

YES 

3.5 

3.5 

2-A 

NO 

NO 

YES 

2/2.5 

2.25 

2-C 

YES* 

NO 

NO 

2/1/1/1 

1.25 

2-1 

NO 

NO 

NO 

1/1/1 

1 

2-2 

YES 

NO 

UNSURE 

2/1 

1.5 

2-3 

YES 

YES 

YES 

3 

3 

2-4 

YES 

YES 

YES 

3/2/1 

2 

2-6 

YES 

NO 

YES 

2.5 

2.5 

2-7 

YES 

YES 

YES 

3 

3/3 

2-9 

YES 

YES 

YES 

3 

3 

2-10 

YES 

YES 

YES 

4 

4 

2-11 

YES 

YES 

YES 

3 

3 

3-C 

NO 

NO 

NO 

1/1.5 

1.25 

3-1 

YES 

NO 

YES 

2/3/2 

2.33 

3-2 

YES 

YES 

YES 

3/3 

3 

3-3 

YES 

YES 

YES 

4/3.5 

3.75 

3-6 

YES 

YES 

YES 

3/3 

3 

3-7 

YES 

YES 

YES 

4 

4 

4-C 

NO 

NO 

UNSURE 

1.5/2 

1.75 

4-1 

NO 

NO 

UNSURE 

1 

1 

4-3 

YES 

YES 

YES 

2/3/3 

2.67 

4-4 

YES 

YES 

YES 

3/3/2 

2.67 

4-6 

NO 

NO 

UNSURE 

2 

2 

4-7 

NO 

NO 

UNSURE 

1 

1 

4-10 

YES 

YES 

YES 

4 

4 

4-11 

YES 

YES 

YES 

4 

4 

5-1 

YES* 

NO 

YES 

3/3 

3 

5-3 

NO 

YES 

YES 

1 / 3/3/2 . 5/1 

2.83 

5-4 

YES 

YES 

YES 

2.5 

2.5 

5-5 

YES 

YES 

YES 

3 

3 

5-6 

YES* 

NO 

YES 

3 

3 

5-7 

YES* 

NO 

YES 

3 

3 

5-11 

YES 

YES 

YES 

3.5 

3.5 

6-1 

YES 

YES 

YES 

4 

4 

6-2 

NO 

NO 

UNSURE 

1 

1 

^Indicates  Type  I  PIO  predicted  using  criterion  only 
^PIOR  =  Pilot  Induced  Oscillation  Rating 
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3.  With  the  exception  of  one  configuration,  all  the 
configurations  predicted  to  be  Type  III  PIO  prone  were  also 
predicted  to  be  Type  I  PIO  prone.  What  implication  does  this 
have  on  the  actual  PIO  seen  in  flight? 

Overall,  the  major  drawback  to  Stoith's  theory  appears  at 
this  point  to  be  the  need  for  a  properly  parameterized  pilot 
model.  The  other  problem  (although  admittedly  of  academic 
interest  only)  is  the  possible  physical  discrepancy  that  may 
arise  by  predicting  both  types  of  PIO. 

anith's  theory  predicts  only  whether  or  not  an  aircraft 
has  a  tendency  to  PIO,  but  says  nothing  about  the  severity  of 
the  PIO.  From  the  above  results,  it  appeared  that  anith's 
magnitude  criterion  might  be  useful  as  a  PIO  rating 
predictor.  The  magnitude  of  the  pilot-felt  normal 
acceleration  to  pitch  rate  at  the  subjectively  predictable 
resonant  frequency  (using  the  C  criterion)  is  plotted  versus 
PIO  rating  is  Figure  14.  Only  those  configurations  which 
have  a  phase  margin  of  less  than  15  degrees  for  the 
pilot-felt  normal  acceleration  to  stick  force  loop  are 
plotted.  Suggested  boundaries  for  PIO  ratings  of  one  through 
four  are  indicated  by  the  dashed  lines  and  are  summarized  on 
the  figure.  The  flight  test  portion  of  this  project 
attempted  to  refine  these  boundaries,  as  will  be  discussed 
later . 
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Bandwidth  Method 


Hoh's  proposed  theory  does  not  directly  predict  PIO 
tendency  or  rating,  but  predicts  levels  of  handling 
qualities.  Hoh's  proposed  bandwidth  criterion  boundaries  for 
handling  qualities  levels  are  shown  in  Figure  15.  Hoh 
originally  used  the  LAHOS  pilot  ratings  to  develop  the 
proposed  Category  C  boundaries  in  the  figure.  The  flight 
test  data  frcm  this  study  were  used  to  refine  the  proposed 
Category  C  boundaries,  as  discussed  later. 

Hoh's  bandwidth  method  was  used  to  see  if  any 
correlation  existed  between  PIO  ratings  frcm  the  LAHOS  data 
and  the  bandwidth  criterion.  Figure  15  also  shows  Hoh's 
bandwidth  criterion  for  the  LAHOS  configurations.  The 
average  PIO  rating  is  noted  next  to  each  configuration.  In 
general,  high  PIO  ratings  appear  to  be  associated  with  low 
bandwidths  and  large  phase  delays.  Proposed  boundaries  for 
PIO  ratings  of  one  through  four  are  indicated  in  Figure  16  by 
the  dashed  lines.  The  flight  test  data  frcm  this  study  were 
used  to  refine  these  proposed  boundaries,  as  discussed  later. 
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Figure  15.  Correlation  Between  Hoh's  Bandwidth  Criterion  and  PIO  Ratings,  LAHOS  Data 


IV.  Flight  Test  Method 


The  flight  test  portion  of  this  project  was  conducted  as 
part  of  a  USAF  Test  Pilot  School  systems  project*  known  as 
HAVE  PIO.  The  test  team  consisted  of  three  project  pilots 
and  two  engineers,  including  the  author.  In  addition,  three 
Calspan  safety/instructor  pilots  were  involved.  The  previous 
flying  experience  of  the  three  project  pilots  is  summarized 
in  Table  12. 


Table  12 

Project  Pilot  Experience 


Pilot 

Aircraft 

Hours 

A 

F-15 

1200 

AT-38 

40 

B 

F-lll 

1050 

AT-38 

30 

T-37 

60 

C 

F-4EJ 

700 

T-33 

200 

Based  on  the  previous  analytical  results,  18  different 
aircraft/ flight  control  system  (FCS)  combinations  were 
selected  for  flight  test  on  the  USAF/Calspan  variable 
stability  NT-33.  These  configurations  included  four  sets  of 
Level  1  approach/landing  short  period  dynamics  combined  with 
14  different  flight  control  systems.  The  PIO  tendencies  and 
the  frequencies  of  the  configurations  were  predicted.  Each 
configuration  was  then  flight  tested,  and  the  actual  PIO 
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tendencies  and  frequencies  were  compared  to  those  predicted. 


Test  Item  Description 

The  test  aircraft,  NT-33A  S/N  51-4120,  is  a  modified, 
two  seat,  jet  trainer  owned  by  theUSAF  Flight  Dynamics 
Laboratory  and  operated  by  Calspan  Corporation,  Buffalo,  New 
York.  The  aircraft  is  capable  of  variable  dynamic  response 
and  control  system  characteristics  (Ref  14).  The  variable 
stability  system  (VSS)  modifies  the  static  and  dynamic 
responses  of  the  basic  NT-33A  by  commanding  control  surface 
positions  through  full  authority  electrohydraulic  servos.  A 
programmable  analog  computer,  associated  aircraft  response 
sensors,  control  surface  servos,  and  an  electrohydraulic 
force-feel  system  provides  the  total  simulation  capability. 
Figures  17  and  18  show  a  block  diagram  of  the  variable 
stability  and  a  sketch  of  the  flight  control  system.  The 
instructor/safety  pilot  varies  the  computer  gains  through 
controls  located  in  the  rear  cockpit,  allowing  changes  in 
airplane  dynamics  and  control  system  characteristics  in 
flight. 
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IN  fUGHT  SIMULATION 


Figure  17.  Variable  Stability  NT-33A  block  Diagram 


Figure  18.  Control  System  Layout 
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Test  Instrumentation  and  Data  Reduction 

The  NT-33  instrumentation  system  consisted  of  the 
following  items: 

1.  An  on-board  Ampex  AR  700  magnetic  tape  recording 
system  with  2.25  hours  recording  capability  was  used  to 
record  28  data  parameters  at  100  samples  per  second,  as 
listed  in  Table  13. 

2.  An  AN/ANH-2  voice  recorder  set  manufactured  by  the 
Pierce  Wire  Recorder  Corporation  provided  45  minutes  of 
recording  time  for  interphone  and  UHF  radio  communications. 

3.  A  Head  Up  Display  (HUD)  video  recorder  was  used  to 
record  the  HUD  field  of  view  and  interphone  communications 
for  all  approaches  and  landings. 

After  each  flight  project  pilots  reviewed  their  HUD 
video/audio  and  inflight  pilot  comment  cards  and  summarized 
their  comments  on  a  mission  summary  sheet.  In  addition, 
individual  project  pilot  Cooper-Harper,  PIO,  and  confidence 
rating  factors  were  determined  and  recorded  for  each 
configuration.  To  maintain  standardization,  these  comments 
and  ratings  were  reviewed  to  ensure  project  pilots  used 
similar  criteria  when  assigning  PIO  and  Cooper-Harper 
ratings.  The  inflight  pilot  comment  card,  Cooper-Harper 
rating  scale,  and  Confidence  ratings  are  shown  in  Figures  19, 
20,  and  21.  The  PIO  rating  scale  was  depicted  in  Figure  13. 

The  data  were  analyzed  using  several  qualitative  and 
quantitative  analysis  techniques.  Project  pilot  comments 
were  used  to  qualitatively  describe  the  aircraft  PIO 
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Table  13 


Test  Instrunentation 
NT-33  Digital  Tape  Parameters 


DIGITAL  CHANNEL 
NUMBER 

RECORDED  VARIABLE 

ENGINEERING  UNITS 

1 

Elevator  deflection 

(measured  at  strut) 

degrees 

2 

Elevator  deflection 

(measured  at  surface) 

degrees 

3 

Aileron  deflection 

(measured  at  strut) 

degrees 

4 

Aileron  deflection 

(measured  at  surface) 

degrees 

5 

Elevator  stick  deflection 

inches 

6 

Lateral  stick  deflection 

inches 

7 

Rudder  pedal  deflection 

inches 

8 

Elevator  stick  force 

pounds 

9 

Lateral  stick  force 

pounds 

10 

Rudder  pedal  force 

pounds 

11 

Event  marker 

N/A 

12 

True  airspeed 

feet/second 

13 

Roll  rate 

degrees/second 

14 

Pitch  rate 

degrees/second 

15 

Yaw  rate 

degrees/second 

16 

Normal  acceleration 

(measured  at  c.g.) 

g’s 

17 

Angle  of  attack 

degrees 

18 

Angle  of  sideslip 

degrees 

19 

Pitch  angle 

degrees 

20 

Roll  angle 

degrees 

21 

Normal  acceleration 
(measured  at  pilof  's 

station) 

g's 

22 

Not  used 

23 

Not  used 

24 

Not  used 

25 

Pitch  error 

degrees 

26 

Roll  error 

degrees 

27 

Not  used 

28 

Not  used 
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Feel  System  Characteristics: 

-  Forces/Displacements? 

-  Pitch  Sensitivity? 

Pitch  Attitude  Control? 

-  Initial  Response? 

-  Final  Response? 

-  Predictability? 

-  Any  special  piloting  techniques/compensation  required? 

-  Tendency  toward  PIO? 

Task  Performance: 

-  Airspeed  Control? 

-  Touchdown  Point  Accuracy? 

-  Sink  Rate  at  Touchdown? 

-  Runway  Alignment? 

-  Level  of  aggressiveness  used  to  control  touchdown  point? 

-  Special  control  techniques  required  in  flare? 

-  If  approach  was  abandoned,  was  it  due  to  poor  handling  qualities 
or  severe  PIO? 

Additional  Factors: 

-  Wind/Turbulence 

-  Lateral-Directional  Characteristics 
Summarize  Evaluation: 

-  Major  problems,  good  features 
Review  Ratings: 

-  PIO  Rating,  Cooper-Harper  Rating,  Confidence  Factor 

Figure  19.  In-flight  Pilot  Comment  Card 
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AOEOUACY  FOB  SELECTED  TASK  OR 
REQUIRED  OPERATION* 


DEMANDS  ON  THE  PILOT 

AIRCRAFT  CHARACTERISTICS  •  M  SELECTED  TASK  OR  REOUtREO  OPERATION* 


Cscaiiani 

Pilot  compensation  nol  e  foctor 

for 

Highly  desirable 

dtured  parformenca 

Good 

Pilot  compentotion  nol  a  factor 

for 

Negligible  dtfioancitt  , 

desired  partarmance 

Fou  -  Soma  mildly 

Mmmol  pM  comptnsoiion  nquirtd 

for 

inpitosam  deficiencies 

da«rad  parlormonca 

Mmor  buf  annoying  ' 

Desired  performonce  requires  moderate 

deficiencies 

pilot  compensoiion 

Moderately  oPpciicnabie 

Adequate  performance  requires 

deficiencies 

considerable  pilot  compensation 

Very  objedionotte  bul 

Adequate  performonce  requires  eitensive 

toierobie  deficiencies 

plot  compensation 

Deficiencies 

require 

improvement 


Mojor  deficiencies 


Major  deficiencies 


Major  deficiencies 


Adequate  performonce  nol  otloinoble  with 
monmum  tolerable  pilot  compensation 
Controllability  not  in  oultlion 


Considerable  pilot  compensation  is  required 
tor  control 

Intense  pilot  comptnsoiion  ii  required  to 
retain  control 


.  \ 
ImprovtfDtnt 

Confrol  null  be  lose  dr  mg  some  poriion  of  required 

10  ) 

mondoiory 

— y  Mo,or  deficiencies  ■  ^ 

J 

\  "Definition  of  required  operation  involve!  designation  of  flight  phose  ond/or  subphoses  with 

Plo*  decu -am  I  accompanying  conditions 


DEFINITIONS  FROM  TN-D-5153 
5ATION  PERFORMANCE 


COMPENSATION 
Th*  measure  of  additional  pilot  effort 
and  attention  required  to  maintain  a 
given  level  of  performance  in  trie  face  of 
deficient  vehicle  charactenstics 

HANDLING  QUALITIES 

Those  Qualities  or  characteristics  of  an 
aircraft  that  govern  the  ease  and  preci¬ 
sion  with  which  a  pilot  is  eCHe  to  perform 
the  tas«s  required  m  suppon  of  an  air¬ 
craft  rpie 

MISSION 

The  composite  o<  pilot-vehicle  functions 
mat  must  be  performed  to  fulfill  opeia- 
nonai  lequuements  May  be  specified  lor 
a  role  complete  flight  flight  phase  o' 
flight  subphase 


The  precision  of  control  with  respect  to 
aircraft  movement  that  a  pilot  is  able  to 
achieve  in  performing  a  last  (Pilot- 
vehicle  performance  is  a  measure  of 
handling  performance  Pilot  perform¬ 
ance  is  a  measure  of  the  menne'  o' 
efficiency  with  which  a  pilot  moves  the 
principe1  controls  in  performing  a  task  ) 

ROLE 

The  function  o'  purpose  the:  delmes  me 
primary  use  0*  an  aircraft 

TASK 

The  actual  work  assigned  a  pilot  to  be 
performed  m  completion  of  or  as  repie 
sentative  of  a  designated  flight  segment 


WORKLOAD 

The  integrated  pnysical  ana  mem  4  eh  or  required 
to  perform  a  specified  piloting  las* 


Figuri'  20.  Cooper-Harper  Pilot  Rating  Scale 
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PILOT  CONFIDENCE  FACTORS 


CLASS  A 

A  pilot  may  assign  a  rating  with  a  relatively 
high  degree  of  confidence,  although  he  may  have 
mild  reservations  because  of  incomplete  or  inade¬ 
quate  simulation  of  motion  cues,  disturbances, 
visual  information,  or  other  factors  affecting 
pilot  workload. 

Supplementary  tasks,  if  needed,  can  be 
adequately  provided  by  the  pilot. 


CLASS  B 

A  pilot  can  assign  a  rating  with  only  a  moderate 
level  of  confidence  because  of  uncertainties 
introduced  by  a  lack  of  representative  environmental 
disturbances  as  well  as  incomplete  or  inadequate 
simulation  of  motion  cues,  disturbances,  visual 
information,  or  other  factors  affecting  pilot 
workload. 

Supplementary  tasks  may  be  desired,  but  are  not 
available. 


CLASS  C 

A  pilot  can  assign  a  rating  with  only  minimum 
confidence  because  considerable  pilot  extrapolation 
is  required  due  to  an  incomplete  task,  thereby 
requiring  considerable  reliance  on  self-imposed  tasks 
and  maneuvers  for  assessment. 

This  may  also  be  aggravated  by  incomplete  or 
very  limited  simulations  of  motion  cues,  disturbances, 
visual  information,  or  other  factors  affecting  pilot 
workload . 


Figure  21.  Pilot  Confidence  Factors  Scale 
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tendencies  and  handling  qualities  during  the  landing  task. 
Pilot  comments,  Cooper-Harper  and  PIO  ratings,  and  strip 
chart  data  of  stick  force,  pitch,  pitch  rate  and  normal 
acceleration  were  used  to  determine  whether  or  not  the 
aircraft  had  a  PIO  tendency  during'  the  approach  and  flare. 

For  data  analysis,  a  PIO  was  defined  as  a  sustained  pitch 
rate  oscillation. 

Test  Configurations 

Landing  longitudinal  PIO  tendencies  and  flying  qualities 
were  evaluated  using  four  pairs  of  short  period  natural 
frequency  and  damping  ratio  combined  with  fourteen  different 
flight  control  system  configurations.  All  short  period 
dynamics  met  MIL-F-8785C  Level  1  boundaries  for  the  landing 
approach  (Category  C).  These  configurations  are  depicted  in 
Figure  22  and  are  also  listed  in  Table  14.  Table  14  also 
gives  the  dynamic  characteristics  of  the  14  different  flight 
control  systems. 


77 


Figure  22.  Flight  Test  Configurations. 


Table  15  shows  the  actual  flight  control  system/aircraft 
dynamics  combinations  used.  The  phugoid  and 
lateral-directional  characteristics  were  held  constant  and 
met  MIL-F-8785C  Level  1  criteria.  These  characteristics  are 
listed  in  Appendix  A. 
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Table  14 


NT-33A  Longitudinal  Dynamics  and  Flight  Control  Systems 


Dynamics 

Configuration 


U) 

n 


sp 


2 

3 

4 

5 

2.4 

4.1 

3.0 

1.7 

0.64 

1.0 

0.74 

0.68 

K(s+t  ) 

First  order  systems:  - 

(s+V 


K 


Second  and  fourth  order  systems:- 

/ 


2  2  2 
(s“+2 s+m  )(s  +2^0)-,  s +w  ) 

ii|  ii  i  * ^  2  n  2 
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Table  15 


Flight  Control  System  and 
Aircraft  Dynamics  Combinations 


The  short  period  dynamics  were  set  by  the  NT-33A 
instructor/safety  pilot  by  adjusting  the  appropriate  variable 
stability  gain  controls  in  the  rear  cockpit.  Predetermined 
flight  control  system  characteristics  were  also  selected  by 
the  rear  seat  pilot.  The  stick  force  per  inch  (gearing)  was 
selected  by  the  first  pilot  to  fly  each  configuration.  From 
that  point  on,  the  gearing  for  each  configuration  remained 
fixed.  The  gearings  for  each  configuration  are  listed  on  the 
pilot  comment  cards  in  Appendix  B. 
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Landing  Task 


A  landing  task  was  defined  to  allow  repeatability  of  the 
results.  Each  pilot  flew  up  to  three  approaches  for  each 
configuration  -  a  straight  in  approach  followed  by  two  offset 
approaches,  one  to  each  side  of  the  runway.  After  making  all 
three  approaches,  the  project  pilots  assigned  both  PIO  and 
Cooper-Harper  ratings  as  a  measure  of  PIO  tendency  and  pilot 
performance  and  workload.  For  the  landing  task,  a  PIO  was 
defined  as  a  sustained  oscillation  which  interfered  with  the 
accomplishment  of  the  task  and  required  the  pilot  to  reduce 
his  gain  or  remove  himself  from  the  loop.  A  PIO  tendency  was 
defined  as  an  undesirable  motion  which  did  not  necessarily 
interfere  with  the  accomlishment  of  the  task. 

The  offset  landing  task  for  this  project  was  a  visual 
approach  with  a  lateral  offset  and  a  correction  to  centerline 
prior  to  touchdown.  Figure  23  depicts  the  runway  landing 
task  parameters.  The  size  of  the  lateral  offset  was 
approximately  150  feet.  Due  to  runway  maintenance,  the  left 
150  feet  of  the  300  foot  wide  Runway  22  at  Edwards  AFB  was 
closed  during  the  test  period.  The  centerline  of  the 
remaining  150  foot  wide  runway  was  used  for  touchdown.  The 
aircraft  was  flown  on  the  desired  glidepath  using  the  ILS 
until  the  beginning  of  the  overrun,  then  the  correction  to 
the  desired  touchdown  point  was  initiated.  The  safety  pilot 
assisted  in  maintaining  a  constant  offset  correction  and 
break  point  among  the  three  project  pilots. 
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Figure  23. 


Runway  Landing  Task  Parameters 


The  touchdown  zone  was  1000  feet  long  starting  at  500 
feet  from  the  threshold  and  extending  to  1500  feet  from  the 
threshold.  The  desired  touchdown  aimpoint  was  1000  feet  from 
the  threshold  and  within  5  feet  of  centerline.  Even  though 
the  ILS  glidepath  intersected  the  runway  at  the  desired 
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touchdown  point,  the  pilots  were  still  required  to  make  a 
large  longitudinal  correction  (push  over)  due  to  the  long 
flare  characteristics  of  the  NT-33A.  Each  landing  was 
treated  as  a  "must  land"  situation,  unless  the 
instructor/safety  pilot  or  project  pilot  determined  that 
safety  of  flight  would  be  compromised  in  an  attempt  to  land. 
Table  16  summarizes  the  evaluation  task  performance  criteria 
used  to  assign  Cooper-Harper  ratings  to  the  visual  landing 
task . 


Table  16 


Task  Performance  Standards 


Desired 
No  PIOs 

Tochdown  within  5  ft  of 
centerline  (main  wheels  on 
centerline) 

Touchdown  aimpoint  ±  250  ft 
Approach  airspeed  +  5  kts 


Adequate 

Touchdown  within  25  ft  of 
centerline  (tip  tank  on 
centerline) 

Touchdown  at  aimpoint  +  500ft 
Approach  airapeed  -5/+10  kts 


An  electronic  step  elevator  input  was  accomplished  on 
downwind  before  reconfiguring  to  allow  post  flight  time 
response  verification  of  the  configuration  dynamics  just 


flown . 
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V.  Flight  Test  Results  and  Analysis 


All  planned  configurations  were  flown  by  at  least  two  of 
the  project  pilots.  PIO  ratings,  Cooper-Harper  pilot 
ratings,  and  pilot  comments  were  collected  for  each 
configuration.  All  pilot  comments  are  summarized  in  Appendix 
B.  Representative  time  history  plots  of  F  ,  0,  q,  and  a„ 

S  Z  p 

for  the  last  30  seconds  of  the  approach  and  landing  for  each 
configuration  are  shown  in  Appendix  C.  Table  17  summarizes 
the  Cooper-Harper  pilot  ratings  and  the  PIO  ratings  for  each 
configuration . 


Table  17 

PIO  and  Cooper-Harper  Ratings,  HAVE  PIO  Data 


Configu¬ 

ration 

#  of 
Flights 

PIO 

Ratings 

Pilot 

Ratings 

2-B 

4 

3/2/2/1 

7/3/3/3 

2-1 

3 

1/1/1 

2/2/3 

2-5 

3 

4/4/5 

10/7/10 

2-7 

3 

4/3/2 

7/4/4 

2-8 

3 

4/4/4 

8/10/8 

3-D 

2 

1/1 

2/2 

3-1 

3 

3/2/2 

5/3/4 

3-3 

3 

3/1/1 

7/2/3 

3-6 

2 

2/2 

5/4 

3-8 

3 

4/3/4 

8/5/8 

3-12 

2 

4/5 

7/9 

3-13 

2 

4/5 

10/10 

4-1 

3 

1/1/1 

3/2/3 

4-2 

3 

1/1/2 

3/3/7 

5-1 

2 

1/1 

2/5 

5-9 

2 

4/4 

7/7 

5-10 

2 

5/5 

10/10 

5-11 

3 

2/4/3 

7/7/5 
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Smith's  Theory 


Tables  18  and  19  summarize  the  Type  I  PIO  tendency 
predicted  for  each  configuration  using  the  v  and  z,  criteria, 
respectively.  Table  20  summarizes  the  Type  III  PIO  tendency 
predictions,  and  Table  21  presents  an  overall  summary  of  the 
predictions  using  Smith's  theory.  All  the  configurations 
predicted  to  be  Type  III  PIO  prone  were  also  predicted  to  be 
Type  I  PIO  prone.  For  Type  I  PIO,  the  only  discrepancy 
between  the  two  criteria  occurred  for  Configuration  2-B. 

This  configuration  used  the  lead/lag  pilot  model  and  was 
predicted  to  be  PIO  prone  using  the  z,  criterion  but  not  the  v 
criterion . 

Table  21  also  presents  the  flight  test  PIO  ratings. 

Using  an  average  PIO  rating  of  2  or  greater  as  a  basis  for 
PIO  tendency,  Smith's  theory  correctly  the  PIO  tendency  (or 
lack  of  PIO  tendency)  for  14  of  the  18  configurations.  Two 
configurations  which  were  predicted  to  be  not  PIO  prone  had 
PIO  tendencies  in  flight,  and  two  of  the  configurations  which 
were  predicted  to  be  PIO  prone  had  no  PIO  tendencies  in 
flight . 

Table  22  compares  smith's  predicted  Type  I  and  Type  III 
PIO  frequencies  to  the  actual  PIO  frequency  from  flight.  The 
actual  PIO  frequency  was  determined  from  strip  chart  data  as 
the  average  PIO  frequency  from  all  approaches  for  a 
particular  configuration.  Overall,  it  is  difficult  to  draw  a 
conclusion  about  whether  a  Type  I  or  Type  III  PIO  was 
experienced  because  the  predicted  frequencies  for  the  two 
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Smith's  Type  I  PIO  Predictions  Using  v  Criterion,  HAVE  PIO 
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Table  20 


Smith's  Type  III  PIO  Predictions,  HAVE  PIO  Data 

.  w  i  w  o 

Configu-  R  c  Predir  - 

ration  (rad/sec)  (rad/sec)  PIO? 


2-B 

11.86 

4.67 

NO 

2-1 

7.07 

4.21 

NO 

2-5 

2.39 

2.99 

YES 

2-7 

4.05 

4.19 

YES 

2-8 

3.66 

4.14 

YES 

3-D 

9.09 

4.91 

NO 

3-1 

11.68 

5.00 

NO 

3-3 

5.36 

4.52 

NO 

3-6 

6.90 

4.99 

NO 

3-8 

5.40 

4.94 

NO 

3-12 

2.27 

3.04 

YES 

3-13 

2.95 

3.75 

YES 

4-1 

8.70 

4.65 

NO 

4-2 

5.33 

4.52 

NO 

5-1 

5.79 

3.65 

NO 

5-9 

2.53 

3.40 

YES 

5-10 

2.14 

2.91 

YES 

5-11 

2.93 

3.65 

YES 

^Frequency  at  j^-axis  crossing 
Crossover  frequency  predicted  using  Smith's  formula 
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Table  21 


Smith's  Overall  PIO  Predictions,  HAVE  PIO  Data 

Configu-  Predict  Predict  PIO  in  Actual,  Average 

ration  Type  I  Type  III  Plight  PIOR  PIOR 


2-B 

YES* 

NO 

2-1 

NO 

NO 

2-5 

YES 

YES 

2-7 

YES 

YES 

2-8 

YES 

YES 

3-D 

NO 

NO 

3-1 

NO 

NO 

3-3 

YES 

NO 

3-6 

NO 

NO 

3-8 

YES 

NO 

3-12 

YES 

YES 

3-13 

YES 

YES 

4-1 

NO 

NO 

4-2 

NO 

NO 

5-1 

NO 

NO 

5-9 

YES 

YES 

5-10 

YES 

YES 

5-11 

YES 

YES 

YES 

3/2/2/1 

2 

NO 

1/1/1 

1 

YES 

4/4/5 

4.33 

YES 

4/3/2 

3 

YES 

4/4/4 

4 

NO 

1/1 

1 

YES 

3/2/2 

2.33 

NO 

3/1/1 

1.68 

YES 

2/2 

2 

YES 

4/3/4 

3.68 

YES 

4/5 

4.5 

YES 

4/5 

4.5 

NO 

1/1/1 

1 

NO 

1/1/2 

1.33 

NO 

1/1 

1 

YES 

4/4 

4 

YES 

5/5 

5 

YES 

2/4/3 

3 

♦Indicates  Type  I  PIO  predicted  using  C  criterion  only 
1PIOR  =  Pilot  Induced  Oscillation  Rating 
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types  are  too  close  together.  In  general,  it  appears  that 
the  flight  test  frequencies  are  closer  to  the  Type  I 
predictions.  Two  of  the  configurations  which  were 
incorrectly  predicted  for  PIO  tendency  (3-1  and  3-6)  have  the 
largest  frequency  errors.  Since  the  predicted  resonant 
frequency  is  dependent  upon  the  parameterization  of  the  servo 
pilot  model,  a  "generic"  pilot  model  as  used  in  this  study 
may  not  be  accurate  enough  for  Smith's  theory. 


Table  22 


Smith's  Predicted  PIO  Frequencies 
and  Actual  Flight  Test  PIO  Frequencies 


Configuration 

Type  I 
Frequency 
(rad/sec) 

Type  III 
Frequency 
( rad/ sec) 

Average  PIO 
Frequency 
( rad/sec) 

2-B 

3.9 

_  _  _ 

00 

• 

2-1 

4.6 

— 

— 

2-5 

2.8 

2.4 

2.8 

2-7 

4.1 

4.1 

4.1 

2-8 

3.8 

3.7 

4.3 

3-D 

6.8 

— 

— 

3-1 

7.9 

— 

10.4 

3-3 

5.0 

— 

— 

3-6 

6.2 

— 

8.4 

3-12 

2.6 

2.3 

2.2 

3-13 

3.2 

3.0 

3.2 

4-1 

5.4 

— 

— 

4-2 

4.9 

— 

— 

5-1 

3.9 

— 

— 

5-9 

2.9 

2.5 

3.5 

5-10 

2.5 

2.1 

2.6 

5-11 

3.3 

2.9 

3.6 

Table  23  compares  the 

average  flight 

test  PIO  rating 

each  configuration  to  the  PIO  rating  predicted  using  the 
magnitude  criterion  as  described  previously.  The  average 
difference  is  0.6  PIO  ratings,  with  a  maximum  difference  of 
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Table  23 


Predicted  PIO  Ratings  from  Magnitude  Criterion 
and  Flight  Test  PIO  Ratings  (PIOR) 


Configu¬ 

ration 

Predicted 

PIOR 

Actual 

PIOR 

Average 

PIOR 

A  ! 
PIOR1 

2-B 

2 

3/2/2/1 

2.00 

0.00 

2-1 

1 

1/1/1 

1.00 

0.00 

2-5 

4 

4/4/5 

4.33 

+0.33 

2-7 

2 

4/3/2 

3.00 

+1.00 

2-8 

2 

4/4/4 

4.00 

+2.00 

3-D 

1 

1/1 

1.00 

0.00 

3-1 

1 

3/2/2 

2.33 

+1.33 

3-3 

2 

3/1/1 

1.68 

-0.32 

3-6 

1 

2/2 

2.00 

+1.00 

3-8 

2 

4/3/4 

3.68 

+1.68 

3-12 

4 

4/5 

4.50 

+0.50 

3-13 

4 

4/5 

4.50 

+0.50 

4-1 

1 

1/1/1 

1.00 

0.00 

4-2 

1 

1/1/2 

1.33 

+0.33 

5-1 

1 

1/1 

1.00 

0.00 

5-9 

4 

4/4 

4.00 

0.00 

5-10 

4 

5/5 

5.00 

+1.00 

5-11 

3 

2/4/3 

3.00 

0.00 

difference  in  magnitude  between  average  flight  test  and 
predicted  PIO  rating 
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2.0  PIO  ratings.  Figure  24  depicts  a  plot  of  average  PIO 
rating  versus  the  predicted  magnitude  criterion.  Like  the 
LAHOS  data,  there  is  a  definite  upward  trend  between  the 
magnitude  criterion  and  flight  test  PIO  rating.  However, 
there  is  too  much  scatter  in  the  data  at  the  lower  PIO 
ratings  to  draw  any  conclusion  about  a  definite  correlation. 
This  scatter  is  due  both  to  the  small  sample  size  (scatter  in 
PIO  ratings)  and  the  impreciseness  of  the  pilot  model 
(predicted  frequency  affects  the  magnitude  criterion) .  A 
magnitude  criterion  of  greater  than  0.022  g/deg/sec 
accurately  predicts  a  PIO  rating  of  4  or  5,  but  more  data  are 
needed  to  further  define  the  lower  PIO  rating  boundaries. 

It  also  should  be  noted  that  one  of  the  configurations 
with  the  largest  difference  between  predicted  and  actual  PIO 
rating  (3-1)  has  a  very  high  predicted  resonant  frequency 
associated  with  a  very  negative  phase  margin.  This  anomaly 
was  discussed  previously  for  one  of  the  LAHOS  configurations 
(4-C).  HAVE  PIO  Configuration  3-1  supports  the  idea  that 
configurations  with  very  high  frequencies  which  fail  the 
magnitude  criterion  may  still  be  PIO  prone. 

An  attempt  was  made  to  further  analyze  the  flight  test 
data  using  a  frequency  response  analysis  to  produce  power 
spectral  densities  (PSDs)  of  pitch,  pitch  rate,  and 
pilot-felt  notmal  acceleration  using  approximately  the  last 
10  seconds  of  flight  prior  to  touchdown  (or  safety  pilot 
assumption  of  controls).  It  was  hoped  that  this  technique 
would  produce  more  accurate  dominant  PIO  frequencies  than 
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Figure  24.  Correlation  Between  Magnitude  Criterion  and 
Flight  Test  PIO  Rating 
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could  be  ascertained  using  strip  chart  data.  A  further 
attempt  was  made  to  compare  the  dominant  frequencies  to  those 
predicted  for  both  Type  I  and  Type  III  PIO  to  determine  which 
type  had  actually  occurred.  However,  most  of  the  PSDs  did 
not  show  any  particular  dominant  frequencies.  Figure  25 
presents  typical  PSDs  for  a  configuration  which  received  a 
PIO  rating  of  3.  The  only  configurations  which  consistently 
produced  dominant  frequencies  using  the  PSD  analysis  were 
those  which  had  a  sustained  PIO;  i.e.,  a  PIO  rating  of  4  or 
5.  Figure  26  presents  typical  PSDs  for  a  configuration  of 
this  type. 

The  PSD  analysis  was  an  attempt  to  gain  some  physical 
insight  into  the  PIO  phenomenon,  but  was  unfortunately 
inconclusive.  A  better  approach  in  the  future  might  be  to 
use  simulator  studies  in  combination  with  in-flight 
simulations.  By  definition,  it  should  be  possible  to 
duplicate  a  Type  III  PIO  in  a  simulator  but  not  a  Type  I  (due 
to  the  lack  of  normal  acceleration  cues) .  By  comparing  PIOs 
experienced  in  the  simulator  with  those  experienced  in  flight 
for  a  given  configuration,  it  may  be  possible  to  ascertain 
which  type  actually  occurred. 
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Figure  25.  Power  Spectral  Densities  of  Selected  Parameters,  Configuration  2-7,  Average  PIO  Rating 


Figure  26.  Power  Spectral  Densities  of  Selected  Parameters,  Configuration  2-5,  Average  PIO  Rating  ■  4.33 


Bandwidth  Method 


Figure  27  presents  Holt's  proposed  handling  qualities 
level  1,  2,  and  3  boundaries.  Also  shown  are  the  test 
configurations  along  with  their  associated  pilot  ratings. 
These  data  are  also  summarized  in  Table  24.  Hoh's  theory 
correctly  predicted  the  level  of  handling  qualities  for  13  of 
the  18  configurations.  Three  of  the  configurations  which 
were  not  correctly  predicted  are  very  close  to  the  Level  2 
boundary  (2-8,  3-8,  and  5-11).  If  the  boundary  were  shifted 
down  slightly,  these  configurations  would  have  been  correctly 
predicted  as  being  Level  3.  The  other  two  configurations 
which  were  not  correctly  predicted  were  2-B  and  3-D.  These 
configurations  were  both  Level  1  in  flight,  but  were 
predicted  to  be  Level  2.  This  suggests  that  Hoh's  dashed 
line  for  the  Level  1  boundary  could  be  drawn  to  the  6  rad/ sec 
bandwidth  point  instead  of  the  5  rad/sec  point.  Based  on 
this  limited  flight  test  data.  Figure  28  presents  suggested 
boundary  changes  for  Hoh's  bandwidth  criterion. 

Proposed  PIO  rating  boundaries  based  on  the  preliminary 
research  presented  earlier  are  depicted  in  Figure  29,  and 
each  flight  test  configuration  is  plotted  to  predict  its  PIO 
rating.  The  figure  also  includes  the  flight  test  PIO 
ratings.  These  data  are  also  summarized  in  Table  25.  The 
average  difference  in  PIO  rating  is  0.5  with  a  maximum 
difference  of  1.3.  Using  this  limited  flight  test  data,  the 
suggested  PIO  rating  boundaries  were  refined  as  shown  in 
Figure  30. 


97 


I 


o 

in 


o 

0 

CM 


O 

0 

o 


(oss)  Aejap  ssEiy 


98 


Bandwidth  (rad/sec) 

Figure  27.  Bandwidth  Method  Predicted  Handling  Qualities  Levels  and  Flight  Test  Pilot  Ratings 


lilt 
no  v  in 

O'  O'  O'  O' 


U  M  IW 

gill 

©D  <\  0 


«/ 


e&SEt 

f 


C^/ 

AS  / 


r-  c 
o»->.  — 

i  p* 

c  m—  r" 
«-  *  v 
rr^fn  in  © 
hO  ^  I  — 


i‘  - 
hr.  ^ 

—  QCv.  ^ 

syjsw 

Sd 

F«  1  o 

rl 

^  \ 

[O  . 


in  m 
<T^»ru  ^ 
IM  n 


\  5 

\ 


(39S)  Aex^P  asi?qd 


99 


Figure  28.  Proposed  Boundary  Changes  for  Hoh's  Bandwidth  Criteria 


Bandwidth  (rad/sec) 

Figure  29.  Proposed  PIO  Rating  Boundaries  Using  LAHOS  Data,  with  HAVE  PIO  Flight  Teat  PIO  Ratings 


The  major  difference  between  the  LAHOS  data  and  the  HAVE 
PIO  data  was  that  there  were  no  LAHOS  configurations  which 
received  a  PIO  rating  of  5,  but  there  were  several  HAVE  PIO 
configurations  which  were  rated  as  5's.  Reviewing  the  LAHOS 
time  histories  and  pilot  comments  cards,  and  discussing  the 
flights  with  the  LAHOS  pilots  revealed  that  several  of  the 
LAHOS  PIOs  were,  in  fact,  divergent  and  should  have  received 
5's.  Although  the  difference  is  of  academic  interest  only 
(both  a  4  and  5  are  unflyable)  the  difference  helps  explain 
disrepancies  between  the  two  data  bases.  Also,  the 
difference  between  parameters  which  predict  a  4  or  a  5  might 
be  helpful  in  ascertaining  whether  or  not  the  PIO  rating 
scale  can  be  considered  linear. 


Table  24 


Bandwidth  Method  Predicted  Handling  Qualities 
and  Flight  Test  Pilot  Ratings 


Configu¬ 

ration 

Predicted 

Level 

Pilot 

Ratings 

Actual 

Level 

2-B 

2 

7/3/ 3/ 3 

1 

2-1 

1 

2/2/3 

1 

2-5 

3 

10/7/10 

3 

2-7 

2 

7/4/4 

2 

2-8 

2 

8/10/8 

3 

3-D 

2 

2/2 

1 

3-1 

2 

5/3/4 

2 

3-3 

1 

7/2/3 

1 

3-6 

2 

5/4 

2 

3-8 

2 

8/5/8 

3 

3-12 

3 

7/9 

3 

3-13 

3 

10/10 

3 

4-1 

1 

3/2/3 

1 

4-2 

1 

3/3/7 

1 

5-1 

2 

2/5 

2 

5-9 

3 

7/7 

3 

5-10 

3 

10/10 

3 

5-11 

2 

7/7/5 

3 

101 


Table  25 


Predicted  PIO  Ratings  from  Bandwidth  Criterion 
and  Flight  Test  PIO  Ratings  (PIOR) 


Configu¬ 

ration 

Predicted 

PIOR 

Actual 

PIOR 

Average 

PIOR 

A 

PIOR 

2-B 

2 

3/2/2/1 

2.00 

0.00 

2-1 

1 

1/1/1 

1.00 

0.00 

2-5 

3 

4/4/5 

4.33 

+1.33 

2-7 

3 

4/3/2 

3.00 

+0.00 

2-8 

3 

4/4/4 

4.00 

+  1.00 

3-D 

2 

1/1 

1.00 

-1.00 

3-1 

2 

3/2/2 

2.33 

+0.33 

3-3 

1 

3/1/1 

1.68 

+0.68 

3-6 

1 

2/2 

2.00 

+1.00 

3-8 

3 

4/3/4 

3.68 

+0.68 

3-12 

4 

4/5 

4.50 

+0.50 

3-13 

4 

4/5 

4.50 

+0.50 

4-1 

1 

1/1/1 

1.00 

0.00 

4-2 

1 

1/1/2 

1.33 

+0.33 

5-1 

1 

1/1 

1.00 

0.00 

5-9 

3 

4/4 

4.00 

+1.00 

5-10 

4 

5/5 

5.00 

+1.00 

5-11 

3 

2/4/3 

3.00 

0.00 

difference  in  magnitude  between  average  flight  test  and 
predicted  PIO  rating 
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VI.  Conclusions  and  Recommendations 


Current  PIO  prediction  techniques  predict  only  PIO 
tendency.  This  research  attempted  to  determine  if  existing 
techniques  could  be  used  to  determine  the  severity  of  the  PIO 
as  well;  i.e.,  predict  PIO  rating. 

Smith's  theory  can  accurately  predict  both  PIO 
tendencies  and  frequencies  with  a  properly  parametreized 
pilot  model.  This  need  for  an  accurate  pilot  model  is  a 
major  drawback  to  Snith's  theory.  There  appears  to  be  a 
correlation  between  Snith's  magnitude  criterion  and  actual 
PIO  ratings.  Further  research  needs  to  be  performed  using 
more  accurate  pilot  models  along  with  existing  data  bases  to 
determine  this  correlation  more  precisely. 

Hoh's  bandwidth  method  appears  premising  as  both  a 
handling  qualities  predictor  and  PIO  rating  predictor.  Hoh's 
method  is  particularly  attractive  because  it  requires  no 
pilot  model.  The  data  collected  in  thi3  study  were  used  to 
suggest  changes  to  Hoh's  proposed  level  boundaries  and  to 
draw  proposed  PIO  rating  boundaries  for  the  Category  C  phase 
of  flight.  Finally,  it  was  pointed  out  several  times  that 
there  is  uncertainty  as  to  whether  or  not  the  PIO  rating 
scale  is  linear  with  respect  to  increasing  PIO  severity. 
Defining  the  relationship  between  PIO  rating  and  seme  easily 
measurable  parameter  (such  as  pitch  or  normal  acceleration) 
would  greatly  aid  PIO  research.  Such  a  relationship  would 
help  in  analyzing  PIO  ratings,  pilot  comments,  and  time 
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history  data  to  determine  consistency  among  pilots.  This 
might  also  help  understand  configurations  with  "clifflike" 
qualities;  i.e.,  sane  pilots  have  no  problems,  but  other 
pilots  behave  just  slightly  more  aggressively  and  have  severe 
problems . 

One  important  question  in  PIO  research  remains 
unanswered,  however.  What  really  causes  PIO?  Are  PIOs 
caused  by  the  switch  fran  tracking  pitch  to  tracking 
acceleration  at  certain  frequencies,  as  Smith  proposes?  The 
accurate  PIO  frequency  predictions  in  this  study  certainly 
seem  to  add  weight  to  his  theory.  Or  are  PIOs  simply  caused 
by  large  time  delays,  which  can  cause  the  pilot  to  become  out 
of  phase  with  his  inputs,  as  Hoh's  method  suggests?  This 
study  showed  that  Hoh's  method  can  predict  PIO  tendencies  as 
well  as  Staith's,  but  with  an  entirely  different  premise.  It 
is  possible  that  the  real  cause  is  a  combination  of  the  two 
theories.  Bandwidth  and  phase  delay  may  in  fact  be  related 
somehow  to  Smith's  "switching"  frequency.  Further  research 
using  new  and  existing  PIO  data  bases  and  simulator  studies 
is  warranted  to  determine  a  real  understanding  of  the  PIO 
phenomenon . 
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APPENDIX  A 


ESTIMATED  STABILITY  DERIVATIVES  AND  TRANSFER  FUNCTIONS 
FOR  THE  LAHOS  AND  HAVE  PIO  CONFIGURATIONS 
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Overall  Aircraft  Configuration 

The  aircraft  was  always  flown  in  the  power  approach 
configuration  (gear  down,  flaps  30  degrees,  speed  brake 
extended) .  The  only  variation  among  approaches  was  the 
approach  airspeed,  which  varied  with  aircraft  weight  (fuel 
remaining).  These  varying  approach  speeds  are  indicated 
below: 


Fuel  Remaining  (Gals) 

150 

250 

350 

450 

550 


Approach  Speed  (KIAS) 

125 

130 

135 

140 

140 


A  nominal  touchdown  speed  of  120  KIAS  (U  =205 

o 

feet/second  and  WQ=25  feet/second)  was  used  for  defining  the 
dynamic  characteristics  of  both  the  LAHOS  and  HAVE  PIO 
configurations.  Other  nominal  aircraft  characteristics  are 
listed  below. 

n  fa~  4.5  g/ radian 

Z 

T02  =1.4  second 

lx  =  6.43  feet  (distance  between  center  of  gravity  and 

pilot's  station) 

•v  * 

Phugoid  characteristics: 

=  0.17  radian/ second 

5  =  0.15 

ph 

T =  12  second 
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Lateral-directional  characteristics : 

*  1.3  radian/ second 

=  0.2 

d 

U/6|d  -  1.5 

r  »  75  second  x  -  0.3  second 

s  s 

Feel  system  characteristics! 


Longitudinal : 

6 

ES  = 

F 

ES 

0.125 

[0.6,26] 

in/ lb 

Lateral : 

5 

0.25 

AS  r 

F 

AS 

[0.7,26] 

in/ lb 

Directional : 

6 

0.017 

RP  = 

frp 

[0.7,26] 

in/ lb 

Actuators:  5625.0 

[0.7,75] 

where  2  2 

[?,w  ]  =  (s  +2C^s+w)  and 
n  n  n 

(s)  *  (s+a) 

For  LAHOS,  the  gearing  ratio  between  the  elevator  and 
the  stick  position  was  selected  by  the  pilot  for  each  flight 
evaluation  of  a  configuration.  For  HAVE  PIO,  the  gearing 
ratio  was  selected  for  each  configuration  by  the  first  pilot 
to  fly  it;  on  subsequent  evaluations  the  gearing  ratio 
remained  set  at  the  value  initially  selected. 
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LAHOS  Configurations 

The  dynamics  characteristics  for  the  LAHOS 
configurations  are  shown  in  Table  26. 


Table  26 


LAHOS  Dynamic  Characteristics 


6-1 


Parameter  1-1 

2-1 

3-1 

4-1 

5-1 

( YF-17 ) 

0) 

sp 

1.03 

2.30 

2.19 

2.00 

3.90 

1.94 

^sp 

0.73 

0.57 

0.25 

1.06 

0.53 

0.65 

X 

u 

-0.041 

-0.041 

-0.041 

-0.041 

-0.041 

-0.041 

Xw 

0.11 

0.11 

0.11 

0.11 

0.11 

0.11 

xq 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X6e 

0.0032 

0.0032 

0.0032 

0.0032 

0.0032 

0.0032 

Zu 

-0.25 

-0.26 

-0.26 

-0.26 

-0.26 

-0.26 

Z 

w 

-0.75 

-0.75 

-0.75 

-0.75 

-0.75 

-0.75 

Z 

q 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Z6e 

1.1 

1.1 

l.l 

1.1 

1.1 

1.1 

M 

u 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

M 

w 

-0.00232 

-0.01875 

-0.02239 

-0.00663 

-0.05934 

-0.01184 

M 

q 

-0.76 

-1.83 

-0.29 

-3.49 

-3.25 

-1.75 

“6e 

0.33685 

0.33685 

0.33685 

0.33685 

0.33685 

0.33685 

©0 

4.5 

4.5 

4.5 

4.5 

4.5 

4.5 

Note  that  there  are  no  values  for  Z,  and  M,  .  These 

w  w 

derivatives  could  not  be  determined  directly  with  the 
parameter  identification  technique  used  by  Calspan.  Z.  was 


109 


assumed  to  be  zero,  and  is  effectively  included  in  the 
derivatives  above.  For  example,  the  Mg  listed  above  is 
really  the  basic  Mg+UQM^  and  Mw  is  the  basic  i.e., 

and  M^  are  considered  “lunped"  stability  derivatives 
(6:211),  as  used  in  Section  II. 

The  transfer  functions  for  each  of  the  configurations 
were  obtained  using  equations  [6],  [7],  [9],  and  [10]  and  are 
summarized  below. 

LAHOS  Configuration  1-1 
A  -  [0.17, 0.13][0. 73, 1.03] 

0 

=  0.33685(0.0827)  (0.7007) 

N^ZP=  -1.066(0) (0.0266) [0.05, 6. 85] 

6  e 

Nw  =  -0.0037(0)[0.23,0.18] 

8 

LAHOS  Configuration  2-1 
A  =  [0.15,0.17][0.57,2.30] 

e 

>I6e  =  0.33685(0.0848)  (0.6950) 

N^ZP=  -1 .066(0) (0.0260) [-0.03, 6.83] 

e0 

N  =  -0.0089(0) (0.0817) (-1.359) 

8 

LAHOS  Configuration  3-1 
A  =  [0.13,0.20][0.25,2.19] 

n!  =  0.33685(0.0850) (0.6929) 

6e 

Njzp=  -1 .066(0) (0.0260) [0.09, 6.82] 

0  e 

N  9  =  -0.0014(0) (0.0389) (-15.1) 

w8 
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LAHOS  Configuration  4-1 

A  =  [0. 25, 0.12](1. 391) (2.841)  «  [0 . 25, 0. 12][1 . 06, 2 .0] 
N®e  =  0.33685(0.0848) (0.6946) 

N?eP“  -1*066(°H0*0260)C“0*15»6*833 

6  e 

Nu  =  -0.0170(0)[0.79,0.21] 

WR 


LAHOS  Configuration  5-1 
A  =  [0.15,0.18][0.53, 3. 89] 

N®  0.33685(0.0859) (0.6813) 

°e 

N^zp=  -1.066(0)  (0.0259) [-0.14, 6. 77] 
e 

N#  *  -0.0159(0)  (O  ^iO..)  (-2.8940) 

o 


LAHOS  Configuration  6-1 

A  =  [0.15,0.16][0.65,1.94] 

9 

N6  =  0. 33685(0. 0853)(0. 6870) 

e 

N?eP=  -1*066(0) (0.0258) [-0.02, 6.80 j 

e 

Nw  =  -0.0085(0) (-0.0051) (-0.6528) 

8 


The  flight  control  systems  for  the  LAHOS  configurations 
are  shown  in  Table  27,  and  Table  28  shows  the  actual  flight 
control  system/aircraft  dynamics  combinations  used. 
Configuration  6-1  was  the  original  YF-17  configuration  and 
6-2  was  the  modified  YF-17  configuration.  These  two 
configurations  are  described  in  the  main  body  of  the  report. 
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Table  27 


LAHOS  Flight  Control  Systems 


A 

B 

FIRST 

C 

ORDER 

1 

FILTERS 

•  2 

3 

4 

5 

K 

2.5 

3.0 

5.0 

1.0 

10.0 

4.0 

2.0 

1.0 

T 

1 

4.0 

3.33 

2.0 

0.0 

— 

— 

— 

— 

T 

0 

10.0 

10.0 

10.0 

0.0 

10.0 

4.0 

2.0 

1.0 

SECOND  AND  FOURTH  ORDER  FILTERS 


6 

7 

8 

9 

10 

11 

K 

256 

144 

81 

36 

16 

65536 

1 

0.7 

0.7 

0.7 

0.7 

0.7 

0.93 

“"1 

16 

12 

9 

6 

4 

16 

C2 

_  —  _ 

—  “  “ 

“  “  — 

” 

™  — 

0.38 

W 

n  2 

—  —  “ 

—  —  — 

“  “ 

16 

First  order  systems;  Kfs+j^) 

<S+V 

Second  and  fourth  order  systems;  K 

(s2+2cu>  s+w  2)(s2+2CU?  s+W  2) 

nt  n  j  n  2  n  2 
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Table  28 


LAHOS  Plight  Control  System  and 
Aircraft  Dynamics  Combinations 


HAVE  PIO  Configurations 

In  order  to  determine  approximate  stability  derivatives 
for  the  HAVE  PIO  configurations,  LAHOS  l-l  was  used  as  a 
baseline  configuration,  and  the  feedback  characteristics  of 
the  aircraft  were  used  to  estimate  the  new  stability 
derivatives.  The  three  stability  derivatives  which  can  be 
modified  are  Zw,  “w  ,  and  Mg  as 


Z'  =  Z  -  Z  §  Ka/U 
w  w  °e  a/  o 

M'  ■  M  -  M<5  Ka/U  [27] 

w  w  e  07  o  L  J 

M '  =  M  -  M  5  K  /U 
q  q  e  q  o 


The  gains  Ka  and 
desired  stablity 


Kg  are  the  feedback  gains  needed  to  get  the 
derivatives.  Using  the  short  period 
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approximation,  the  stability  derivatives  in  turn  determine 
the  short  period  damping  ratio  as 


a)2  =  z 1 M 1  -  M'U 

n  w  q  wo 

2?U)n  =  -IT  -  Z^  [28] 


Substituting  LAHOS  1-1  values  in  [27]  and  then 
substituting  the  result  into  [28]  gives: 


small 
* 


of  =  1.04587  +  0.25264K  +  0.34093K*  +  0.0Q*MK„K 

n  q  “  /  “  q  £29] 


1.51  +  0.33685K  +  l.lKa/205.0 

q 


[29]  can  then  be  solved  for  and  Ka  for  a  given  value 
of  x,  and  u)n.  Kq  and  Ka  are  then  substituted  back  into  [27] 
along  with  LAHOS  1-1  values  to  give 


=  -0.75  -l.lKa/205.0 

■  -0.0023213  -  0. 3368 5KQ/ 205.0  [30] 

M'  =  -0.76  -  0.33685K 

q  q 

Equation  [30]  provides  estimates  of  the  stability 

derivatives  required  to  give  the  desired  values  of  w  and  S. 

n 

Note  that  Ka  and  are  meaningless  except  as  a  stepping 
stone  to  determine  the  new  derivatives. 
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The  above  technique  was  tested  using  the  remaining  five 
LAHOS  configurations.  Using  LAHOS  1-1  as  a  baseline 
configuration,  the  technique  was  able  to  estimate  the 
stability  derivatives  for  the  remaining  LAHOS  configurations 
within  five  percent. 

The  technique  was  then  applied  to  the  desired  short 
period  characteristics  of  the  HAVE  PIO  configurations.  The 
dynamic  characteristics  of  the  HAVE  PIO  configurations  are 
shown  in  Table  29.  Note  that  HAVE  PIO  configuration  2-1  is 
very  close  to  LAHOS  configuration  2-1,  and  HAVE  PIO 
configuration  5-1  is  very  close  to  LAHOS  configuration  6-1 
( YF-17 ) . 
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Table  29 


HAVE  PIO  Dynamic  Characteristics 


Parameter 

2-1 

3-1 

4-1 

5-1 

w 

sp 

2.41 

4.22 

3.04 

1.70 

^sp 

0.63 

0.97 

0.73 

0.68 

X 

u 

-0.041 

-0.041 

-0.041 

-0.041 

X 

w 

0.11 

0.11 

0.11 

0.11 

X 

q 

0.0 

0.0 

0.0 

0.0 

X5e 

0.0032 

0.0032 

0.0032 

0.0032 

Zu 

-0.26 

-0.26 

-0.26 

-0.26 

Zw 

-0.80642 

-0.92116 

-0.84168 

-0.76979 

Z 

q 

0.0 

0.0 

0.0 

0.0 

z  & 

°e 

1.1 

1.1 

1.1 

1 . 1 

M 

u 

0.0 

0.0 

0.0 

0.0 

M 

w 

-0.01960 

-0.05474 

-0.03040 

-0.00838 

M 

q 

-2.26560 

-7.27889 

-3.59834 

-1.54220 

M6e 

0.33685 

0.33685 

0.33685 

0.33685 

00 

4.5 

4.5 

4.5 

4.5 
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The  transfer  functions  for  each  of  the  configurations 
were  obtained  using  equations  [6],  [7],  [9],  and  [10]  and  are 
summarized  below. 

HAVE  PIO  Configuration  2-1 
A  =  [0.15,0.17][0.63,2.41] 

0 

N  *  =  0.33685(0.0845)(0.6990) 

u  e 

N§ZP=  -1.063(0) (0.026) [-0.06, 6.86] 

0 

Nwg  -0.0111(0) (0.0108) (-1.019) 


HAVE  PTO  Configuration  3-1 
A  =  [0.17,0.16][0.97,4.22] 

0 

N  ^  e  =  0.33685(0.08470(0.6987) 

N|zp=  -1.0626(0) (0.0262) [-0.44, 6. 85] 

N9  =  -0.0355(0) (-0.6566) (-0.0048) 
W8 


HAVE  PIO  Configuration  4-1 
A  =  [0.16,0.16][0.73,3.04] 

6 

N*  =  0.33685(0.0846) (0.6988) 

NjzP=  -1.0626(0) (0.0261) [-0.16, 6.86] 

e 

N9  =  -0.0176(0) (0.0084) (-0.9395) 

W8 


HAVE  PIO  Configuration  5-1 

A  =  [0.16,0. 15][0. 68, 1.70] 

N®e  =  0.33685(0.0845X0.6989) 
e 

N^zp=  -1 .0626(0) (0.0260) [-0.01, 6.86] 

°e 

N9  =  -0.0075(0) (-0.0422) (-0.3432) 

8 
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APPENDIX  B 


PILOT  COMMENT  CARDS 
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This  appendix  contains  the  pilot  comment  cards  for  each 
configuration.  The  cards  include  all  summarized  comments, 
PIO  and  Cooper-Harper  ratings,  and  information  on  the  flight 
control  system  (FCS)  used  and  the  gearing  selected.  The  FCS 
information  is  summarized  in  the  same  manner  as  the  rest  of 
the  report : 

(a)  =  (s+a) 
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119 


s  a 

I  ! 

-  <  Q» 


O  Q.  — 

O 

6 


K> 

K> 

3  s 

N  6  Ki 


t  a.  a> 
u 
4iA  u. 


i  0 

i 

<N 

i  * 

So 
S- 
5° 
5  £ 

o  3  x 
— ■  ii  a 

c  8i 

u.  U  13 


0) 


32 

C  2  o  f 

ss!h 


si 


ft! 

2  2  g  • 

yil 

T3  T9  m  Q 

2235: 


Z 

O 


!  fc  z  > 

J  JUU 

ss&sig 

S**fc°5 

ftij-So 

M  *  *.  A  A  A 


o 

in 


o  • 

-*  in  — 
m  4«  m 


Ifllll 

</>  •  <*  5  ?  o  a 

532225= 


0^0 
rt  3  r; 


< 

at 


-Sa-5u’*e  * 

ySos-SSj,. 

S&5*5£dS 

is  (V  ”>  sr  T>  iu  < 
®io2UOo,u  A 
S'  «r  i  ■  <v  i  -  <  55  a;  5 

“■•••••■a  p- 


ui 

CJ  ui 


cc  tv* 

Z 

I 


n  o 

ui  i 


ii 

1 1 

-  <-> 

^  « 

&  3 

1  a 

ii 

o 

22 

2  3 

W  o 

c  w 


h  e 

§2 

°  a. 

o  31 

5  a 


££ 


T’CVI  — 


Z  .. 

a  (A 

|§ 

is 

5  at 

go. 

ui  a 
IM  O 


T  • 

r 


n 

o 

o 

CO 


120 


> 


> 


£  « 

£  u< 

< 

N 

1  5 

s 

R  <  2 

DATE: 

PILOT. 

IP 

at 

UJ 

£ 

5 

u 

w  5 

O  — 

»•  M  It 

w 

w»sp 

?5P 

FCS 

II 

CD 

z 

H 

< 

CK 

z 

CJ 

1  NO.: 

2-1 

5  M 

g20 

*- 

*“ 

O  1-  ° 

"  j  II 

H 

CD 

z 

o  1  2 

*  i  § 

5 

at 

tj  O  tu 
u.  a  o 

o 

a. 

8c 


«* 

f 

| 

o 

m  n  in 

M  M  « 

4= 

Ci 

£ 

O 

is 

ill 

o  ©  o 

«*  •  «f 

L.  k.  4. 

i 

D  € 

22 

Iliii 

Ill 

11 

|l|S8g! 

s|liSSi 


E  j  *  al  5  9 

,*;nt 


•  ao'oiij 

H  T  w-i  _ 

•  CJ  </> 

I  < 


*  S“  X  <  I  i 

J  T  7  ^  O  ...  , 

®  —  O  =>  O  O  0.  ui 

jjj  <  I  •  (»•  I--  <  in  o: 


a  a 
U  o 


2 

QC 

O 

14. 

uj  a 
cj  uj 
z  o. 

^  (K 
^  o 


Z 

2  * 


z  .. 

O  01 

££ 

1b 

I  “ 


<Y  |-  — 


UJ  Of 

r\j  o 


|j? 


r 

r 


122 


•3GDD  FEATURES.  Good  flying  qualities  eGGOO  FEATURES:  Quick  response.  Corrections  were  quickly  damped 

Felt  like  normal  (tore  except  for  quicker  response 
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tc  low  frequency,  medium  amplitude  PIO.  PlO's  up  and  Medium  frequency,  low  amplitude  motion,  not  divergent 

away  on  downwind 

•GOOD  FEATURES: 
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SOX  FEATURES:  »3G0D  FEATURES:  Aircraft  flaw  wall  both  final  end  flare 

Aircraft  responds  wall  to  pilot 
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•3033  FEATURES:  No  lag  or  delay  was  noted 

The  motion  was  fast  with  no  residual  overshoots 
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•3C3D  FEATURES:  Overall.  no  delay  or  lag  was  noted.  «G000  FEATURES:  Good  flying  qualities. 

Pilot  felt  in  phase  with  input 
Intentional  inputs  during  flare  did  not  induce  PIO. 
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FIT  NO  =  17  RUN  NO:  4  M*p  :  3  0  OAT£:  29MAYB6  FLTK3  :  4  RUN  NO.:  I  co*p  :  3.0  DATE:  l7nAV06 

CONFIGURATION  N0.»  4-1  %SP  *  0.74  PILOT:  B  CONFIGURATION  NO.*  4-1  .  0.74  PILOT.  C 

GEARING:  0.200  FCS  *  I  IP:  HARPER  GEARING:  0.200  FCS  «  I  IP:  EASTER 


5*5 

..  _I  c 

M 

Cl 

£ 

a  t 

E  > 

1  i 

^  <n<* 

2  o  5 
Egt 


2  225. 

•O  «  —  —  c 

S 

iiiiL 

sellEB 

3  -j  a  h  u  z 

cl'^°g 

*  ^  iT  ft.  a.  0. 

*••••• 


s§i 

■*  m  "  _ 
m  «<  m  a 
« ~  ♦.  & 
w  •  w 

Ills  I 

sisifis 

dd 

UOZ-lSjn 

«oci<S5u? 

7“o>qS 

«n*l<xu5^5 

2 "  *  y  2  y  s  a  $ 


|ll 

e  |! 
e  £  • 
..  <  3  o 


W  _ J  — *  I  I  M  3 

y  &  a  !?  far*  V  Si  j 

as&a  Iv.-m  s? 

•  •  •  •  h  •  •  •  X  m 

—  v 


O  TJ  O  C 

SS3SS 


— •  0  a* 

^  S  W 

m  m  m  £ 

♦I  N  «  ft 

W  OB 

ft  o  o  c 

«  a  n  5 

t  U  L  2  V 

iJi  5  51  J  c  c  < 

0  0  @  Q  OS 

a  o  o  ”  E  z  z 


o 

Oft*" 
—  O  o 
(S  C  Z 


01  * 
I  }*{  I 

h  u  u 

Ege 

>•  u  O. 


£ui  i-«S 
a  i  t  V  . 

IS&Sgi 


□  _j  h  M  2 

-  1  <  5  6  “ 
£  £°2 
«  -  ».  a.  n.  o. 

as  •  •  •  •  • 


o  5  § 

-gSS 

u  p  z  = 

5  u  ^  ^ 

«  o  o  ^ 
TuiQ> 

g£x  < 

2  £  y  £ 

*  -  a  ^ 

<v 

•  •  • 

* 


2  ! 


•  w  a  il 

*>S^ 

qJJ 

I  Ui  2  O 
LI  Q£  u  Ol 
n  io  u  < 
u  po.  ui 
i  <  w  <v 


|  y  a?  /- 

i  >  ,J  •*. 

N-  •  •  • 


o  o> 

r  c  o) 

5  u  uj 


«  r  <n 

r  •  • 


137 


3CX  FEATURES.  Good  control  harmony  end  coordination 
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Aircreft  initial  response  tended  to  be  slow. 
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•  30 3D  FEATURES  •GOOD  FEATURES:  PIO  quickly  Hits  out  with  rtductd  inputs 
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APPENDIX  C 


SELECTED  PARAMETER  TIME  HISTORIES 
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TIME  -  SEC 

Figure  31.  Selected  Parameter  Time  Histories,  Configuration  2-B 


Figure  32. 
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TIME  -  SEC 

Figure  34.  Selected  Parameter  Time  Histories,  Configuration 


Figure  37.  Selected  Parameter  Time  Histories,  Configuration 


Figure  38.  Selected  Parameter  Time  Histories,  Configuration 
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Figure  39.  Selected  Parameter  Time  Histories,  Configuration 
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Figure  43.  Selected  Parameter  Time  Bistories,  Configuration 
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Figure  44.  Selected  Parameter  Time  Histories,  Configuration  4-2 


NT-33  ICNGITUDINAL  PIO  PREDICTION  TEST  SHORT  PERIOD  FREQUENCY  =  1.7  RAD /SEC 

CONFIQJRATICN  5-1,  PILOT  A,  FLIGHT  3  SHORT  PERIOD  DAMPING  RATIO  =  0.68 

POWER  APPROACH  CONFIGURATION,  130  KIAS  ECS  =  NONE 
12,900  LB  GROSS  WT,  C.G.=25.3%  MAC 


Figure  45.  Selected  Parameter  Time  Histories,  Configuration 


NT- 33  LONGITUDINAL  PIO  PREDICTION  TEST  SHORT  PERIOD  FREQUENCY  =  1.7  RAD/SEC 

CONFIGURATION  5-9,  PILOT  C,  FLIGHT  5  SHORT  PERIOD  DAMPING  RATIO  =  0.68 

-20.04-  POWER  APPROACH  CONFIGURATION,  135  KIAS  PCS  =  36 . 0/ [S2-I-2 ( . 7)  (6)s+62] 

13,550  LB  GROSS  WT,  C.G.=25.7%  MAC 
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Figure  47.  Selected  Parameter  Time  Histories,  Configuration  5-10 
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